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C^l ' We describe a semiclassical method to calculate universal transport properties of chaotic cavities. 

, While the energy-averaged conductance turns out governed by pairs of entrance-to-exit trajectories, 

the conductance variance, shot noise and other related quantities require trajectory quadruplets; 
' simple diagrammatic rules allow to find the contributions of these pairs and quadruplets. Both pure 

symmetry classes and the crossover due to an external magnetic field are considered. 

^ ■ PACS numbers: 73.23.-b, 72.20.My, 72.15.Rn, 05.45.Mt, 03.65.Sq 

'"ctf ' I- INTRODUCTION 

Mesoscopic cavities show universal transport properties - such as conductance, conductance fluctuations, or shot 
noise - provided the classical dynamics inside the cavity is fully chaotic. Here chaos may be due to either implanted 
impurities or bumpy boundaries. A phenomenological description of these universal features is available through 
, random-matrix theory (RMT) by averaging over ensembles of systems (whose Hamiltonians are represented by ma- 
trices) pj. For systems with impurities, one can alternatively average over different disorder potentials. However, 
experiments show that even individual cavities show universal behavior faithful to these averages. 

In the present paper we want to show why this is the case. To do so we propose a semiclassical explanation 
of universal transport through individual chaotic cavities, based on the interfering contributions of close classical 
trajectories. This approach generalizes earlier work in 0, 0, 0, and is inspired by recent progress for universal 
q ' spectral statistics [g, 0, lS Oi- Our semiclassical procedure often turns out to be technically easier than RMT; 
transport properties are evaluated through very simple diagrammatic rules. 

We consider a two-dimensional cavity accommodating chaotic classical motion. Two (or more) straight leads are 
attached to the cavity and carry currents. We shall mostly consider electronic currents, but most of the following 
. ideas apply to transport of light or sound as well, minor modifications apart. 
\Q ' The leads support wave modes ("channels") e lkXi cos9i sin(kyi sin \0i\); the subscripts i = 1,2 refer to the ingoing 
and outgoing lead, respectively, Xi and y^ with < yi < Wi are coordinates along and transversal to the lead. Here, 
uii is the width of the lead, k the wave number, and 9i the angle enclosed between the wave vector and the direction 
of the lead. Dirichlct boundary conditions inside the lead impose the restriction kwi sin = a^7r with the channel 
index <2j running from 1 to TV,;, the largest integer below Classically, the <Zj-th channel can be associated with 

trajectories inside the lead that enclose an angle 0i with the lead direction, regardless of their location in configuration 
space. The sign of the enclosed angle changes after each reflection at the boundaries of the cavity, and angles of both 
' signs are associated to the same channel. 

We shall determine, e.g., the mean and the variance of the conductance as power series in the inverse of the number 
of channels N = Ni + In contrast to much of the previous literature, we will go to all orders in -4. We shall be 
interested both in dynamics with time reversal invariance ( "orthogonal case" ) and without that symmetry ( "unitary 
O ■ case" ) . For electronic motion time reversal invariance may be broken by an external magnetic field. For that latter 
case, we shall also interpolate between both pure symmetry classes by account for a weak magnetic field producing 
• magnetic actions of the order of U. 
k> | We will always work in the semiclassical limit, and thus require the linear dimension L of the cavity to be large 
compared to the (Fermi) wavelength A. When taking the limit j- — > 0, the number of channels N cx j- (w ~ ~ W2) 
will be increased only slowly. The width of the openings thus becomes small compared to L. For this particular 
semiclassical limit, the dwell time of trajectories inside the cavity, Tjj oc ^ grows faster than the so-called Ehrenfest 
time Te oc In k- oc In Ti. Interesting effects arising for T E /T D of order unity [Hill El 

are thus discarded. 

Following Landauer and Biittiker [Til Il5l ] , we view transport as scattering between leads and deal with amplitudes 
taia 2 f° r transitions between channels ai and a 2 . These amplitudes form an N± x N 2 matrix t = {t aia2 }- Each t aia2 
can be approximated semiclassically, by the van Vleck approximation for the propagator, as a sum over trajectories 
connecting the channels a\ and a 2 |l6|. 

t aia2 ~-±= J2 A a c^' n . (1) 
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The channels exactly determine the absolute values of the initial and final angles of incidence 0\ , 62 of the contributing 
trajectories; again both positive and negative angles are possible. In QJ, T H denotes the so-called Heisenberg time, 
i.e., the quantum time scale l-nfip associated to the mean level density /?■ The Heisenberg time diverges in the 
semiclassical limit like Tjj = 2nhp ~ (2^ft)J ' -i with the volume of the energy shell and / the number of degrees of 
freedom; we shall mostly consider / — 2. The "stability amplitude" A a (which includes the so-called Maslov index) 
can be found in Richter's review [Tg. Finally, the phase in Q depends on the classical action S a = J p ■ dq. 
Within the framework just delineated, we will evaluate, for individual fully chaotic cavities, 

• the mean conductance (ti(tt^)} (Actually, the conductance is given by -%Ytr(#t)\, taking into account two 

2 

possible spin orientations; we prefer to express the result in units of ■%), 

• the conductance variance ((ti tt^) 2 ) — (tr(tt^)\ , 

• the mean shot noise power 

(tr(fcfct _ ttktf)), in units of ^- i J 

• for a cavity with three leads, correlations between the currents flowing from lead 1 to lead 2 and from lead 1 to 
lead 3 depending on the corresponding transition matrices >2 ), i' 1 ^ 3 ) as (tr(t( 1— >2 )f( 1_>2 ) f( 1_>3 )i( 1_>3 ) 

• the conductance covariance at two different energies which characterizes the so-called Ericson fluctuations. 

Here the angular brackets signify an average over an energy interval sufficient to smooth out the fluctuations of the 
respective physical property. We will see that after such an averaging each of these quantities takes a universal form 
in agreement with random-matrix theory, without any need for an ensemble average. To show this, we shall express 
the transition amplitudes as sums over trajectories as in Q). The above observables then turn into averaged sums 
over pairs or quadruplets of trajectories, which will be evaluated according to simple and universal diagrammatic 
rules. We shall first derive and exploit these rules for the orthogonal and unitary cases and then generalize to the 
interpolating case (weak magnetic field). 

Due to the unitarity of the time evolution, we could equivalently express transport properties through reflection 
amplitudes and trajectories starting and ending at the same lead. For the average conductance we have checked 
explicitly that the same result is obtained, meaning that our approach preserves unitarity. 



II. MEAN CONDUCTANCE 



We first consider the mean conductance and propose to show that individual chaotic systems are faithful to the 
random-matrix prediction 0, IT(| 

(G(E)) = (tr(tS)) = {^ L uniter y case (2 ) 
w " \ y >l orthogonal case . v ' 

In the semiclassical approximation ifTJl. the average conductance becomes a double sum over trajectories a, (3 connecting 
the same channels a± and 0,2, 



\ai,a2 a,/3:ai — >a2 



Due to the phase factor e I ( s °~' s ' 3 )/' i , the contributions of most trajectory pairs oscillate rapidly in the limit % — * 0, 
and vanish after averaging over the energy Systematic contributions can only arise from pairs with action differences 
AS = S a — Sp of the order of h~. 



A. Diagonal contribution 



The simplest such pairs involve identical trajectories a. = (3, with a vanishing action difference 0, ITt| . These 
"diagonal" pairs contribute 

<^)>| dia =^/e E i^A ( 4 ) 
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FIG. 1: Scheme of a Richter/Sieber pair. The trajectories a (full line) and /3 (dashed line) connect the same channels a\ and 
02, and differ only inside a 2-encounter (in the box). A Poincare section V intersects the encounter stretches at the points x-pi 
and xp2 in phase space whose configuration-space location is highlighted by two dots. V divides the encounter in two parts 
with durations t s ~ j In ^ and t u ~ jlno. The relevant trajectories are in reality much longer than depicted here; in the 
absence of a potential they consist of a huge number of straight segments reflected at the boundary. 



The foregoing single-trajectory sum may be evaluated using the following rule established by Richter and Sieber 0: 
Summation over trajectories connecting fixed channels is equivalent to integration over the dwell time T as 

]T K| 2 = d Te~^ T = ^L. (5) 



a:ai — >d2 



N rp 

Here, the integrand e t h can be understood as the survival probability, i.e., the probability for the trajectory to 
remain inside the cavity up to the time T. The factor ^ = 2 p( Wi + w ' 2 ) [ s the classical escape rate. Due to fl oc L 2 , 

that rate is proportional to ^ if p is scaled according to p oc L; inversion yields the typical dwell time To — ^S- oc — 
mentioned in the introduction. 

Finally summing over all N\ possible choices for a\ and over the N2 possibilities for 02, one finds [3.ll7| 

(M« t )>U = ^- (6) 

Eq. © reproduces the RMT result for the unitary case, and gives the leading term in the orthogonal case. 



B. Richter/Sieber pairs 



For the orthogonal case, Richter and Sieber attributed the next-to-leading order to another family of trajectory 
pairs. In the following, we shall describe these pairs in a language adapted to an extension to higher orders in j^. 
In each Richter/Sieber pair (see Fig.^l, the trajectory a contains a "2-encounter" wherein two stretches are almost 
mutually time-reversed; in configuration space it looks like either a small-angle self-crossing or a narrow avoided 
crossing. We demand that these two stretches come sufficiently close such that their motion is mutually linearizable. 
Along a, the two stretches are separated from each other and from the leads by three "links" 1 . The partner trajectory 



1 In our previous papers 9, 25, 38] we used the term "loop" to refer to the comparatively long orbit pieces connecting encounter stretches 
to one another or to the openings. We decided to replace this expression by "link" which is more appropriate since the beginning and 
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f3 is distinguished from a only by differently connecting these links inside the 2-encounter. Along the links, however, 
(3 is practically indistinguishable from a; in particular, the entrance and exit angles of a and [3 (defined by the in- and 
out-channels) coincide 2 . The initial and final links are traversed in the same sense of motion by a and f3, while for 
the middle link the velocities are opposite. Obviously, such Richter/Sieber pairs a, (3 can exist only in time-reversal 
invariant systems. The two trajectories in a Richter/Sieber pair indeed have nearly the same action, with the action 
difference originating mostly from the encounter region. 

We stress that inside a Richter/Sieber pair, the encounter stretches and the leads must be separated by links of 
positive durations £1,^2, £3 > 0. For the "inner" loop with duration ti the reasons were worked out in previous 
publications dealing with periodic orbits and [H H(j for more complicated encounters): Essentially, (i is 

obtained from a by switching connections between four points where the encounter stretches begin and end; to have 
four such points the stretches must be separated by a non- vanishing link. The fact that the duration of the initial and 
final links is non-negative (the encounter does not "stick out" through any of the openings) is trivial in the case of 
the Richter/Sieber pair: Since the encounter stretches are almost antiparallel a trajectory with an encounter "sticking 
out" would enter and exits the cavity through the same opening and thus be irrelevant for the conductance. 

Encounters have an important effect on the survival probability 0. The trajectory a is exposed to the "danger" 
of getting lost from the cavity only during the three links and on the first stretch of the encounter. If the first 
stretch remains inside the cavity, the second stretch, being close to the first one (up to time reversal) must remain 
inside as well. If we denote the duration of one encounter stretch by t cnc , the total "exposure time" is thus given 
by T exp = ti + 1 2 + £3 + t cnc ; it is shorter than the dwell time T which includes a second summand i enc representing 

N rp N rp 

the second encounter stretch. Consequently, the survival probability e t h oxp exceeds the naive estimate c t h .In 
brief, encounters hinder the loss of a trajectory to the leads. 

To describe the phase-space geometry of a 2-encounter, we consider a Poincare section V orthogonal to the first 
encounter stretch in an arbitrary phase-space point x-pi. This section must also intersect the second stretch in a 
phase-space point xp>2 almost time-reversed with respect to x-pi. In Fig. ^ the configuration-space locations of 
x-pi and x-p2 are highlighted by two dots. For a hyperbolic, quasi two-dimensional 3 s yste m, the small phase space 
separation between the time-reversed Txp 2 of x-p 2 and x-pi can be decomposed as [2lll22ll23| 

Tx P2 - x-pi = se s (x-pi) + !/e"(x P i) , (7) 

where e s (x-pi) and e u (x-pi) are the so-called stable and unstable directions at Xpi. If V moves along the trajectory, 
following the time evolution of xpi, the unstable component u will grow exponentially while the stable component s 
shrinks exponentially. For times large compared with the ballistic time (L/v with v the velocity) the rate of growth 
(or shrinking) is given by the Lyapunov exponent A (not to be confused with the wavelength also denoted by A), 

u(t) ~ w(0)e At 

s(t) ~ s(0)e- At . (8) 

By our definition of a 2-encounter, the stable and unstable components are confined to ranges — c < s < c, 
— c < u < c, with c a small phase-space separation. The exact value of c will be irrelevant, except that the transverse 
size of the encounter in configuration space, ~ c/VmX with m the mass, must be small compared with the opening 
diameters. It should also be small enough to allow mutual linearization of motion along the encounter stretches. As a 
consequence, the time between V and the end of the encounter is t u ~ j In t^t, i.e., the time the unstable component 
needs to grow from u to ±c. Likewise the time between the beginning of the encounter and V reads t s ~ j In A . 
Both times sum up to the encounter duration 

1 c 2 

tenc = tu + ts ~ V In I 1 ■ (9) 

A \su\ 

A glance at Fig.^shows that the times t-p\,t-p2 of the piercing points xpi, xp2 (measured from the beginning of the 
trajectory) are now given by 

t-pi = tl +t u , tj>2 =tx + t C nc + t2 +t s . (10) 



the end of such a piece may be far removed from each other (in the case of the initial and final link and of links between different 
encounters). 

2 Following Richter and Sieber we find the semiclassical estimate for a conductance component {|ta 1 a 2 | 2 } between t wo g iven in- and 
out-channels and demand therefore that all contributing trajectories have the same in- and out-angles. An alternative |12| is to replace 
summation over channels in the formula for the transport property by integration; then the channel numbers of the contributing 
trajectories found through an additional saddle point approximation will not be i nteg er. 

3 Our treatment can easily be extended to / > 2, see Il9l and the Appendices of Ist l20l . 



5 



Finally, the stable and unstable coordinates determine the action difference as [22, (see also 0,^3) 

AS = su. (11) 
The encounters relevant for the transport phenomena have action differences of order ft and thus durations t enc ~ 

-i 2 -i 2 

j In T^gj of the order of the Ehrenfest time Tg = j In 

With this input, we can determine the average number of 2-encounters inside trajectories a of a given dwell time 
T. In ergodic systems, the probability for a trajectory to pierce through a fixed Poincare section V in a time interval 
{t-p2, tp2 + dt-p2) with stable and unstable separations from x-pi inside (s, s + ds) x (u, u + is uniform, and given 
by the Liouville measure -^dt^dsdu. To count all 2-encounters inside a, we have to integrate this density over tj>2 
(to get all piercings x-p 2 through a given "P) and over t-pi (to get all possible x-pi and thus all possible sections P). 
When integrating over tpi, we weigh the contribution of each encounter with the corresponding duration i cnc , since 
the section V may be placed at any point within the encounter; therefore we must subsequently divide by i C nc- The 
integration over the piercing times tpi, t-p 2 may be replaced by integration over the link durations ii,t 2 , which as 
we stressed, must be positive; in addition t% = T — ti — t 2 — 2i cnc must also be positive. Altogether, we obtain the 
following density of stable and unstable coordinates, 



w(s,u)= dhdh—— r. (12) 

A.+^Sw Vt enc (s,u) 

This density is normalized such that integration over all s, u belonging to a given interval of AS* = su yields the 
number of 2-encounters of a giving rise to action differences within that interval. 

To find what Richter/Sieber pairs contribute to the conductance (j2J|, we replace the sum over (3 by a sum over 
2-encounters inside a or, equivalently, an integral over w(s, u). The additional approximation 4 Ap w A a yields 

( tr ( ttt ))l R s-^( E j dsdu E \A a \ 2 w(s,u)e isu A . (13) 

\ai,ci2 a:ai^a2 I 

Next, we employ the Richter/Sieber rule to do the sum over a by integrating over the dwell time T, with the integrand 
involving the (modified) survival probability c T a cxp = c T H yi 2 3 cnc; . The integral over T may be transformed 
into an integral over the duration of the final link. Moreover, summation over all channels a% = 1 . . . N\, a 2 = 1 • • • -/V 2 
yields a factor N1N2. We are thus led to 

<tr(«t))| = ^1 / r dtldt2 dt 3 f dsdu ai \ e-^t'i-^-H.+u.(.,.)) e u»/A _ (14) 
This integral factors into three independent integrals over the link durations, 

dUe th = —^j- , (15) 
and an integral over the stable and unstable separations within the encounter, 

1=1 f dsdu- - t Ws ' u) e ls "/ R \ . (16) 

\J nt cnc { s ,u) I v ' 

The encounter integral / can be evaluated if we expand the exponential as e Tj^tocis.u) _ ^ _ ^.^^ «) + .... As 

shown in 0, the constant term yields (dsdu Qt 1 ( SM ) e ' s "^) = (^TT^li) wn i cn oscillates rapidly as ft — > and 
therefore vanishes after averaging. In the semiclassical limit, the value of / is solely determined by the linear term for 
which the denominator t onc (s,u) cancels out, 

I= -^{l dsduelsu,n )=-w (17) 



4 For long trajectories, the derivative {J^~^j m A a is proportional to the so-called stretching factor A Q , i.e., the factor by which 

an initial separation along the unstable direction grows until the end of the trajectory a. This factor can be written as a product of the 
(time-reversal invariant) stretching factors of the individual links and encount er stretches . Since j3 contains practically the same links 
and stretches, we have Ap Si A a . All other factors almost coincide as well (see 18, 20, 23] for the Maslov index), entailing Ag ~ A a - 
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FIG. 2: Families of trajectory pairs (ct,(3) differing in one 2-encounter (a), two 2-encounters (b-f) or in one 3-encounter (g-j). 
In contrast to Fig. 0the cavity is not depicted, and the initial and final points of the trajectories are joined together (junction 
symbolized by two dots and intervening bar, -|-). One orbit pair may result from joining beginning and end of different trajectory 
pairs (like cd, ef, and hij). Arrows indicate the directions of motion inside the encounters, and highlight those links which 
are traversed by a and (3 with opposite sense of motion. Note that all families apart from b) and g) involve almost mutually 
time-reversed encounter stretches and thus require time-reversal invariance (TRI). Together with the diagonal pairs, the ones 
shown here reproduce the mean conductance up to order N1N2/N 3 . 

where we use J dsdue lsu / h — > 2nh and T# = -^j-; all further terms vanish compared to the linear one like A ij;™ c ~ ^ 
and may thus be neglected, given our previous definition of the semiclassical limit. Since all occurrences of Th in Eqs. 
(|14fl . I|15|) and l|17|) mutually cancel, we can formulate the following "diagrammatic rule": Each link yields a factor 
jr and each encounter a factor — N . The result still has to be multiplied with the number of channel combinations 
NiN 2 . Altogether, the contribution of Richter/Sieber pairs to the average conductance is hence determined as 



(tr(ttt))| RS = 



N1N2 
N 2 



(18) 



Our present treatment differs from Richter's and Sieber's original paper 0] in two points: First, one has to exclude 
encounters which stick out of the opening. (In 0, encounters were described through self-crossings in configuration 
space, and only that crossing, approximately in the center of the encounter, had to be located inside the cavity.) 
Second, one must take into account that encounters hinder the escape into the leads. While these two corrections mu- 
tually cancel for Richter/Sieber pairs they will presently turn out of crucial importance for higher-order contributions 
to the mean conductance, as well as for other obscrvablcs like shot noise. 



C. Diagrammatic rules for all orders in 



To proceed to all orders in jj, we must consider pairs of trajectories differing by their connections inside arbitrarily 
many encounters. Each of these encounters may involve arbitrarily many stretches. We shall speak of an /-encounter 



FIG. 3: Reconnections inside a 3-encounter leading to one trajectory (dashed) and one periodic orbit (dotted), rather than a 
single connected partner trajectory. 



whenever I stretches of a trajectory come close in phase space. In time-reversal invariant systems we must also allow 
for encounters whose stretches are almost mutually time-reversed. As a consequence, the resulting conductance will 
depend on the symmetry class: We shall see all higher-order contributions to mutually cancel in the unitary case but 
to yield the 1/N expansion of the RMT result (J2J) in the orthogonal case. 

A list of examples is displayed in Fig. [21 where for later convenience we did not draw the cavity and formally joined 
the initial and final points of the trajectories together. These examples illustrate the simplest among infinitely many 
topologically different families of trajectory pairs. 

The individual families are characterized (i) by the numbers vi of l-encounters in which the two partners differ. 
These numbers can be assembled into a "vector"-? = (v2, «3, W4, . . .), and determine the overall number of encounters 
V(v) = X)i>2 v i an( i the- total number of encounter stretches L(v) — ^2i>2^ v l- The num ber of links exceeds the 
number of stretches by one and reads L(v) + 1. Further characteristics of our families of trajectory pairs are (ii) the 
order in which the encounters are traversed by the trajectory a, (hi) the mutual orientation of the encounter stretches 
(i.e., =E vs. ii, or 3 vs. 3), and (iv) the reconnections leading to the partner trajectory (3. We stress that f3 
must be a single connected trajectory; reconnections leading to, e.g., one trajectory and one periodic orbit as in Fig. 
13 must be excluded. 

All families of trajectory pairs contribute to the conductance according to the same rules as do Richter / Sieber 
pairs: Each link yields a factor -h, each encounter gives a factor —N, and we have to multiply with the number of 
channel combinations N\N2- To prove this assertion, we place a Poincare section V a (with = 1 . . . V) across each of 
the V encounters. Similar as for Richter /Sieber pairs, we characterize each ^-encounter by I — 1 stable coordinates s G j 
(with cr = 1 . . . V and j — 1 ... I — 1), and by I — 1 unstable coordinates u a j 0, measuring the phase-space separations 
between the points where the I encounter stretches pierce through V a - 

All V encounters are thus characterized by J2i>2^ — ^) Vl = L — V stable coordinates, and the same number of 
unstable coordinates. As shown in |9(, these coordinates determine the action difference as AS" = Y) c j s a jU a j. Each 
encounter lasts as long as the absolute values of all coordinates remain below the bound c. Consequently, the duration 
of an encounter is determined by the largest stable and unstable coordinates, the first to reach c. In analogy to 10 , 
the cr-th encounter thus has the duration 

1 c 2 

iencM~-ln : — : ■ (19) 

A maxj |s CT j| x max,,-/ lito-j'l 

Again, the trajectory a may get lost from the cavity only during the links and on the first stretch of each encounter. 
If it survives on that first stretch, it cannot escape on the remaining stretches of the same encounter, since these 
are close to the first. The exposure time T exp is thus obtained as the sum of all link and encounter durations 

?oxp = ^2f=i U + 2cr=i Cc- The survival probability e t h Toxp results, again larger than the naive estimate e t h T . 

We proceed to investigating the statistics of encounters for a given family of trajectory pairs. Generalizing the 
treatment of Subsection III Bl we first keep all V Poincare sections V a fixed: each V a is placed orthogonal to a at a 
phase-space point traversed at a fixed time. The further L — V piercings through these sections may be considered 
statistically independent; the probability density for their occurrence at given times and with given stable and unstable 
coordinates thus reads ^r= v ■ To account for all encounters, we must integrate ql-v over all times for the L — V 
later piercings. Moreover, to include all possible V a , we must integrate over the times of the V phase-space points 
chosen as the origin of a section V a . Since all sets of piercings within the duration t^ nc belong to the same encounter, 
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the latter integral weighs each encounter with a factor t° nc - Exactly as for Richter/Sieber pairs, this factor must 
subsequently be divided out. To simplify our calculation, wc replace the integral over altogether L times of piercings 
by an integral over the durations ti, t^, ■ ■ ■ , tt, of all links except the final one. This replacement is permissible because 
all piercing times may be written as functions of t±, ti, ■ ■ ■ , £l, with the Jacobian of the transformation equal to unity. 
Here, the duration of the final link does not show up, since that link does not precede any encounter stretch 
or piercing point. The integral goes over positive ti, by the same reasoning as for 2-encounters 5 ; moreover, the 
cumulative duration of the first L links and all encounter stretches must be smaller than the dwell time T, to allow 
for a non-vanishing final link. We thus obtain the following density of stable and unstable coordinates 

w(s,u)=[ dh...dt L ^ . (20) 

J ti>0 n L - v Y\ v t a (s u) 

T. L 1 t,+Y. v !*» <t " 11(7=1 encv ' / 



The weight w(s,u) is normalized similarly as in Subsection III Bl Integration over Suj,Uaj corresponding to a given 
interval of AS = Y) c j s a jU a j leads to the number of partner trajectories (3 of a given a, with action difference inside 
that interval, and with the pair (a, (3) belonging to the family considered. 

We can now evaluate the contribution of an arbitrary family of trajectory pairs to the average conductance @ • We 
again approximate Ap « A a , and write the sum over f3 as an integral over w(s, u), 

( tr (tt t ))| fam -^(E Jd L - v sd L - v u \A a \ 2 wMe^^ s ^A . (21) 

As before, the sum over a\,a2 leads to a multiplication with the number of channel combinations N1N2. The sum 
over a can be done using the (modified) Richter/Sieber sum rule, and leads to integration over the duration tL+i 

of the final link, with an integrand involving the survival probability e t h oxp = e T H w^t=i » 2^=1 -/^ Together 
with the integrals over the remaining L link durations in u>(s, u), Eq. (|20|l . all links are now treated equally, and we 
obtain 



is. II) 



i . x NlN2 L +l ( r _ N .\V/ f eiE, wj e -^*^ (a 1 
(tr(ttt))| fem = Yl i^j dUe t h u J Y[\J ds al ...ds a ^ 1 du al ...du a ^ 1 ^ ^ ^ I . (22) 

Just like in (1141) the integral factorizes into several integrals, one for each link and each encounter: Each link gives 

and the encounter integral is determined by the linear term in the series expansion of the exponential e T n * onc< - s,tl ' . The 
encounter integral is slightly changed because the piercing probability A , and the simple integral J dsdue lsu ^ h — > 2nh 

now both appear in the (I — l)-fold power. Each encounter thus yields — ^ (^fp)' 1 = — ^r- Again, all powers of 

Th mutually cancel. We thus find the same diagrammatic rules as above with a factor jj from each of the L + l 
links, a factor — N from each of the V encounters, and a factor N1N2 representing the possible channel combinations. 
(We note that these rules have a nice analogy to our previous work on spectral statistics, see Appendix [DJ . The 
contribution of each family is therefore given by 

(M« t )>l^ = (-i) v( ^^|^ T r- (23) 

To obtain the overall conductance, we must sum over all families. If we let Af(v) denote the number of families 
associated to v, Eqs. © and imply 

W«'» = ^ (1 + U-^^vm^i) - ^ (1 + 1 ■ (") 

Here, each coefficient c rn is determined by the families with given m = L — V , 



L ( v ) — V ( v) — m 

(-l) v( ^Af{v). (25) 



E 



5 Note that the initial and final links must have positive duration also in case they lead to stretches of a parallel encounter. Otherwise 
all parallel stretches of that encounter would enter or leave the cavity through the same lead. This is impossible, since the trajectories 
cannot enter or leave the cavity several times. 
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Our task has thus been reduced to counting families of trajectory pairs and evaluating c m . 

For the lowest orders m, the counting is easily done. We have already seen that for time-reversal invariant systems 
the next-to-leading contribution to conductance originates from the Richter/Sieber family of trajectory pairs differing 
in one 2-encounter, see Fig.Q]or[21i. This family has L = 2, V = 1 thus m = 1. Hence, it gives rise to a coefficient 
Ci = — 1 in the orthogonal case. In the unitary case the absence of such trajectory pairs implies c\ = 0. 

The following coefficient C2 is determined by pairs of trajectories which either differ in two 2-encounters (i.e., L = 4, 
V = 2) and thus contribute with a positive sign, or differ in one 3-encounter (i.e., L = 3, V = 1) and contribute 
with a negative sign. The relevant families are sketched in Fig. [21 In the unitary case, there is only one family 
of the first type (Fig. |5Jd), and one family of the second type (Fig. [2k). Both contributions mutually cancel, i.e., 
c 2 = 0. In the orthogonal case, we must allow for encounters with mutually time-reversed stretches. We then find five 
families contributing with a positive sign (Fig. |2t>-f ) and four families contributing with a negative sign (Fig. E^-j). 
All contributions sum up to c 2 = 1. 

The higher coefficients c m require more involved combinatorial methods to which we now turn. 



D. Combinatorics 

In H, we found a systematic way for summing contributions of families of trajectory pairs that differ in arbitrarily 
many encounters. In particular, we obtained a recursion for the numbers Af(v). Since the treatment of [9j was geared 
towards spectral statistics of closed systems, it was formulated for pairs of periodic orbits rather than pairs of open 
trajectories. It can, however, be easily adapted to open trajectories. We just have to turn trajectory pairs into orbit 
pairs by joining the initial and final points as in Fig. [21 or cut orbits in order to form trajectories. 

The topology of orbit pairs (A, B) was described by "structures ". To define these structures, we numbered the 
encounter stretches of A in their order of traversal, starting from an arbitrary reference stretch; the links of A were 
numbered as well, with the first link preceding the first encounter stretch. Then, each structure is characterized by 
(i) a vector v as above, (ii) a way of distributing the numbered stretches among the encounters, (iii) fixing the mutual 
orientation of stretches inside each encounter, and (iv) shifting connections to form a partner orbit B. 

With this definition, each family of trajectory pairs corresponds to one structure of orbit pairs. We only have 
to glue together the initial and final points of the trajectories as in Fig. [21 and keep the first stretch of a (the one 
following the initial point) as the first stretch of A. Thus, each of the pictures in Fig. [21 represents one structure of 
orbit pairs, and the numbers of orbit pair structures and of trajectory pair families both equal Af(v). 

To illustrate this relation, we consider the two families of trajectory pairs depicted in Fig. [2J? and f. If we join the 
initial and final points for any of these families, we obtain orbit pairs of one and the same topology. But still the 
resulting structures are different, because in Fig. [2^ the first encounter stretch (the one following the initial point of 
the trajectory) precedes a stretch of the same encounter, i.e., the two antiparallel encounters respectively involve the 
stretches (1, 2) and (3, 4). A different choice of the initial stretch as in Fig. [2f means that the two encounters comprise 
the stretches (1,4) and (2,3). Indeed, structures of orbit pairs and families of trajectory pairs are one-to-one. 

For later convenience, we refer to the structure involving one antiparallel encounter (Fig. [2k- ) as the Sieber/Richter 
structure (it was proposed by these authors in see also [24|), to the structure involving two parallel 2-encounters 
(Fig. [21)) as ppi, to the two structures involving a parallel and an antiparallel 2-encounter (Fig. [2t and d) as api, and 
to the structures involving two antiparallel encounters (Fig. [21? and f) as aas [25(. The structure involving one parallel 
3-encounter (Fig.|2fe) will be called pc, and the three structures of the type i± (Fig. Gh-j) ac. 

To formulate our recursion for Af(v), we now denote by Af(v, I) the number of structures of orbit pairs (or families 
of trajectory pairs) related to v, for which the first stretch belongs to an l-encounter. We had established the identity 

A/K0 = ^A/-(fl) (26) 
L(v) 

which has the following intuitive interpretation (not to be confused with the proof in 9J): The probability that the 
first stretch forms part of an ^-encounter is given by the overall number Ivi of stretches belonging to ^-encounters, 
divided by the overall number L{v) of stretches in all encounters; to obtain Af(v,l), we have to multiply Af(v) with 
that probability. Incidentally, the definition ofAf(v,l) implies 5Zi>2-^(^0 = N(v). 
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The relevant recursion relation reads Q 6 

N(v 1 2)-Y J N{v^ k+1 \k + l)= fl-iWtfPM), (27) 

k>2 ' 

with (3 = 2 and ,9 = 1 respectively referring to the unitary and orthogonal case. The symbol yi k > 2 ^ k + 1 \ denotes the 
vector obtained from v if we reduce and v 2 by one, and increase Vk+i by one; likewise v^ 2 ^ is obtained from v if 
we reduce v 2 by one. In general, the list on the left-hand side of the arrow contains the sizes of "removed" encounters, 
whereas the right-hand side contains the sizes of "added" encounters. 

To turn (|27() into a recursion for the coefficients c m , we multiply with (— 1)^W and sum over all v with fixed 
to = M(v) = L{v) - V(v), 

M(v)=m M(v)=m , . M{v)=m 



E (-l) v ^Af(v,2)-J2 E (-l) V W(#' a ^ fc+11 ,fc+l)= (|-l) E (28) 

t? fc>2 ^ ' t? 

In each of the foregoing sums we can replace the summation variable v by v' = yi k ' 2 ^ k + i \ or v" = v^ 2 ^ . Given 
that trajectory pairs associated with v' have one encounter and one encounter stretch less than those of v, we then 
have to sum over v' with M(v') — L(v') — V(v') — (L(v) — 1) — (V(v) — 1) = to. By definition, we should restrict 
ourselves to v' with v' k+1 > 0; however, that restriction may be dropped since v 1 with v' k+1 — have Af(v' , k+ 1) = 0. 
In contrast, trajectory pairs associated to v' have one encounter and two stretches less than those associated to 
v. The pertinent sum runs over v" with M{v") = L(v") — V(v") = (L(v) — 2) — (V(v) — 1) = m — 1. Using 
= _(_i)V(«7') = _(_i)V(^') we can rewrite as 

M(v)=m M(v')=m . . A/(t?")=m-l 

E (-i) l/(i7) ^2) + E E (-i)^'M^fc+i) = -(^-i) E (-ir^Wo- (29) 

fc>2 v' VP / 0// 

The left-hand side now boils down to v^^ =ro (-l)^) ^ fc > 1 7V(t/, fc + 1) = £^ 00=m (-l) v ( fl W(u) = c m while the 
right-hand side reads — ^ — lj c m _i. We thus end up with a recursion for c m , m > 2, 

2 \ [O unitary case ,„„. 

p / I— c m _i orthogonal case . 

For the unitary case we conclude that all off-diagonal contributions to the average conductance mutually cancel; 
the remaining diagonal term, Nl ^ 2 reproduces the random-matrix result. For the orthogonal case an initial condition 
is provided by the coefficient c\ = —1, originating from Richter/Sieber pairs; hence c m = (—1)™. The anticipated 
mean conductance @ is recovered through l|24|) as the geometric series 

^1^2 A , ^ (-l)" l \ ^1^2 



AT \ ^— ' iV TO / N + 1 

\ m— 1 / 

We have thus shown for both symmetry classes that the energy-averaged conductance of individual chaotic cavities 
takes the universal form predicted by random-matrix theory as an ensemble average. 

III. CONDUCTANCE VARIANCE 

Experiments with chaotic cavities also reveal universal conductance fluctuations. In particular, the conductance 
variance agrees with the random-matrix prediction Q 

SggW+i) uni ^ case (32) 

N(N+ir(N+3) orthogonal case . 



Eq. 1271 is a special case of Eqs. (42) and (54) in with 1 = 2. To understand the equivalence, note that in Jj we allowed for 
"vectors" v' including a non- vanishing component , which may formally be interpreted as a number of "1-cncounters" . We moreover 
showed that Afffi, 1) = AT(i/[ 1— *i) (see theparagraph preceding Eq. (58) of |3). Applying this relation to v' = {?t 2— '-"J, one sees that 
the A*^ - !?' 2 ^ 1 ' , 1) appearing in Eq. (54) of coincides with Af(v^ 2 ^). 
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Once more, the semiclassical limit offers itself for an explanation of such universality. With the van Vleck approxima- 
tion for the transition amplitudes fl). the mean squared conductance turns into a sum over quadruplets of trajectories, 

<(tr(ttt)) 2 ) = /j2 WUW^) = ( E E A a A*^A* s e^- s ^- s ^A . (33) 

Here a\,c\ — 1 . . . N± and ai^c-i = 1 . . . N2 are channel indices. The trajectories a and /3 lead from the same ingoing 
channel a\ to the same outgoing channel ai 1 whereas 7 and 8 connect the ingoing channel c\ to the outgoing channel 
C2- We can expect systematic contributions to the quadruple sum over trajectories only from quadruplets with action 
differences AS = S a — Sp + 5 7 — Ss of the order of Ti. 



A. Diagonal contributions 

The leading contribution to l|33|l originates from "diagonal" quadruplets with pairwise coinciding trajectories either 
as (a = (3, 7 = 8), or as (a = 8, = 7); both scenarios imply vanishing action differences. The first scenario a = j3, 
7 = 5 obviously leads to f3 connecting the same channels as a, and 8 connecting the same channels as 7, as required 
in H33fl : this holds regardless of the channel indices ffli, Ci, 02, C2. The second scenario (a = 8, j3 = 7) brings about 
admissible quadruplets only if all trajectories connect the same channels, i.e., both the ingoing channels a± — c\ and 
the outgoing channels a-i — ci coincide. The contribution of these diagonal quadruplets to (|33l) may thus be written 
as the following double sum over trajectories a and 7 

((trWlU^ll E KIW+E E KIWV (34) 



= 1 *=f}, ai ->a 2 
c 2 -y = 5:c 1 ^c 2 



The sum over channels just yields the number of possible channel combinations as a factor, namely NfN^ for the first 
scenario and N1N2 for the second one. Doing the sums over a and 7 with the Richter/Sieber rule |J5J we get 



The larger one of the two summands, ^ 2 2 > is cancelled by the squared diagonal contribution to the mean conductance. 



N N 

^■ 2 2 , is cancelled by 1 

In recent works on Ehrenfest-time corrections [T^ . IT^ (which are vanishingly small in our limit Tg <C To) the 
diagonal approximation was extended to include trajectories which slightly differ close to the openings. The relation 
of the methods used in these papers to our present approach is not fully settled yet; further investigation about this 
relation is desirable. 



B. Trajectory quadruplets differing in encounters 

Off-diagonal contributions arise from quadruplets of trajectories differing in encounters; see Fig. 0] for examples. 
Each trajectory pair (a, 7) typically contains a huge number of encounters, where stretches of a and/or 7 come close 
to each other (up to time reversal). Partner trajectories /3, S can be obtained by switching connections within some 
encounters. Together, (5 and 8 go through the same links as a and 7, and traverse each ^-encounter exactly I times, 
just like the pair (a, 7). Consequently, the cumulative action of (/3, 8) is close to the one of (a, 7), with a small action 
difference AS = (S a + S-y) — (Sp + Ss) originating from the intra-encounter reconnections. 

Different quadruplet families are distinguished by the number of Z-encounters, the mutual orientation of encounter 
stretches, their distribution among a and 7, and the reconnections leading to (3 and 8. Similar as the diagonal 
quadruplets, some families of quadruplets involve one partner trajectory whose initial and final links practically 
coincide with those of a and the other one whose initial and final links coincide with those of 7. When these families 
are depicted schematically with encounters suppressed they all look the same and in fact like diagonal quadruplets 
(see Fig.QJi), for which reason we shall refer to them as "d- families" ; examples are depicted in Figs.0J>h. In analogy 
to the diagonal quadruplets, d-quadruplets contribute with altogether N~l iVf + N\N 2 channel combinations. Of these, 
N1N2 arise when a and j3 start and end alike since then the four channel indices involved are unrestricted; when a 
and 8 start and end alike, N1N2 combinations arise since the channels are restricted as a\ = C\, a-i = ci- 

A second type of quadruplet families is drawn schematically in Fig. here one partner trajectory practically 
coincides at its beginning with a and at its end with 7; the other trajectory coincides at its beginning with 7 and 
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FIG. 4: (a) Schematic graph of d- quadruplets, with the hatched area a "black box" containing any number of encounters; one 
of the dashed partner trajectories shares initial and final points with a; the second partner trajectory similarly related to 7. 
(b)-(h) d-quadruplets responsible for the leading-order contribution to the conductance variance. The diagrams (b), (f)-(h) 
containing encounters with antiparallel stretches exist only in the orthogonal case. A diagram may have a "twin" obtained by 
reflection in a horizontal line; the number of symmetric versions of each diagram is indicated by a multiplier underneath, 
(i) Schematic graph of x -quadruplets, encounters suppressed: one of the partners shares initial link with 7 and final link with 
a, the second one connects initial link of a with final link of 7. (j) An ^-quadruplet involving one 2-encounter. 



at its end with a. The simplest example of such an "x-family" involves just one 2-encounter, see Fig. 0,0; not 
surprisingly, our schematic sketch strongly resembles that picture. Since quadruplets contribute to the conductance 
variance only if they connect channels as a, (3 : a\ — > 0,2 and 7, S : c\ — ► C2, ^-families arise only if either the ingoing 
channels or the outgoing channels coincide. If the ingoing channels coincide, a\ = c\, the trajectory coinciding initially 
with 7 and finally with a has the form a\ = C\ — ► 0,2 and may be chosen as (3; the trajectory coinciding initially with 
a and finally with 7 is of the form a\ = c\ — > ci and may be chosen as 5. If 02 = C2, similar arguments hold, with (3 
and 8 interchanged. Thus, x-families arise for iViiVf channel combinations with a\ = c\, and for NfN2 combinations 
with d2 = C2, altogether for NN1N2 possibilities. 

We shall presently find that quadruplet families contribute to the conductance variance according to the same rules 
as do pairs to the mean conductance: Each link yields a factor each encounter a factor — N; moreover, we have to 
multiply with the number of channel combinations, i.e. NfN^ + N1N2 for d-families and NN1N2 for x-families. 

To justify these rules we consider a family with numbers of Z-encounters given by v — (1)2,1)3,1)4,, . . .). Again, v 
determines the total number of encounters V(v) and the number of encounter stretches L(v). The overall number of 
links is now given by L(v) + 2, since there is one link preceding each of the L{v) encounter stretches, and the two 
final links of a and 7 which do not precede any encounter stretch. Similarly as for trajectory pairs, we can determine 
a density w(s, u) of stable and unstable separations; this density will be normalized such that integration over all s, u 
belonging to an interval (AS, AS + dAS) of action differences AS yields the number of pairs j3, S differing from given 
cv,7 such that the quadruplet (a,/?, 7, 5) belongs to a given family and the action difference is inside that interval. 
Using the same arguments as in Subsection III CI one finds w(s, u) as an integral over {Q L ~ V Yl a=1 *cnc( s 7 wl th 
the integration running over the durations of all links, except the final links of a and 7. The integration range must 
be restricted such that all links (including the final ones) have positive durations. To evaluate the contribution of one 
family to the quadruple sum in l|33(l . we may now replace the summation over f3 and 5 by integration over w(s, u), 

<(Mtt t )) 2 >| faro = 2t/E Jd L - V sd L - v u \A a \ 2 \A,\M^u)e^^^^y (36) 

a 2 ,c 2 -y:ci— >c 2 

The sums over a, 7 can be performed using the Richter/Sieber rule to ultimately get further integrals over the 
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durations of the final links of a and 7, with an integrand involving the survival probability exp{— ■^■(^2f=i U + 

Sct=i Cc( s i u ))}- We thus meet with link and encounter integrals of the same type as for the mean conductance. 
All powers of Tjj mutually again cancel, and we are left with a factor from each of the L(v) + 2 links and a factor 

— N from each of the V(v) encounters which altogether give n l(v)-v(v)+2 ■ The summation over a\, c\, a%, c 2 yields the 
number of channel combinations mentioned. 

If we denote by NAv), M x (v) the numbers of d- and i-families associated to v, the sum over all families with fixed 
L(v) — V(v) = m involves the subsums 

L{v)-V{v)=m 
dm = V {-l) V ^N d (v) 



V 

L(v)-V(v)=m 

E 



(~l) v ^ r W x (v) (37) 



which allow to write the yield of all families as 

<(tr(ttt)) 2 ) = { NM + N lN2 ) ( ^ + E ^2 E • (38) 

m— 1 / m—1 



with D arising from <i-familics (including the diagonal contribution) and X from x-families. 

The coefficients d m ,x m are obtained by counting families of quadruplets. That counting is an elementary task for 
small to. The coefficient di accounts for families of <i-quadruplets differing in one 2-encounter (L = 2, V = 1, TO = 1). 
In the unitary case there are no such families, i.e., d\ — 0. In the orthogonal case, we must consider quadruplets with 
one partner trajectory differing from a in a 2-encounter, and one partner trajectory identical to 7, see Fig. 0J>; the 
quadruplet thus contains one Richter/Sieber pair and one diagonal pair. A similar family of quadruplets involves one 
partner trajectory identical to a, and one partner trajectory differing from 7 in a 2-encounter. We thus have d\ = —2. 

The following coefficient di is determined by d-quadruplets differing in two 2-encounters or in one 3-encounter, the 
latter quadruplets contributing with a negative sign. Some of these quadruplets fall into two pairs contributing to the 
average conductance. Quadruplets consisting of one diagonal pair and one pair contributing to the coefficient C2 of 
the average conductance (see Fig- EJ) yield a contribution 2c 2 to d 2 (i.e., in the unitary case and 2 in the orthogonal 
case); the factor 2 arises because either a or 7 may belong to the diagonal pair. In the orthogonal case, there is one 
further family of quadruplets consisting of two Richter/Sieber pairs. Finally, we must reckon with quadruplets that 
do not fall into two pairs contributing to the mean conductance, as depicted in Figs. |3|>e, for the unitary case. Two 
further families are obtained by "reflection", i.e., interchanging a and 7 in Figs. 0Ji andQJj. Taking into account the 
negative sign for Fig.0^ and its reflected version, the respective contributions sum up to 1. In the orthogonal case, 
the additional families in Figs.0f-h and the reflected versions of Fig.Qk and h yield a further summand 1. Altogether, 
we thus find d 2 = 1 in the unitary case and d 2 = 2+1 + 1 + 1 = 5 in the orthogonal case. 

The most important family of ^-quadruplets, see Fig.0], involves a parallel encounter between one stretch of a and 
one stretch of 7. This family, discovered in [3j for quantum graphs, gives rise to a coefficient xi = — 1 for systems 
with or without time-reversal invariance. 

With the coefficients d\,d 2 ,xi, the conductance variance can be evaluated up to corrections of order 
The result 7 

< W « W -<««.))'- M +0 <*) U " i '"' V 7 (38) 

[ ^4 2 + C Ijt ) orthogonal case , 

coincides with the random-matrix prediction (|32|l . We note that Eq. (|39l) could ultimately be attributed only to 
the quadruplets shown in Fig. |IJ>h, since all other contributions mutually cancel. (In particular, the contributions 
proportional to NfN 2 from all (i-quadruplets that consist of two pairs contributing to the conductance are cancelled 
by the squared average conductance. The term proportional to NiN 2 in the diagonal approximation is compensated 
by the contribution of x-quadruplets as in Fig. 



7 In counting orders we assume that all numbers of channels are of the same order of magnitude. 
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To go beyond Eq. (|39|) (and to show that no terms were missed in Eq. 1|39|0 . we must systematically count families 
of d- and ^-quadruplets with arbitrarily many encounters. Similar to the case of conductance this can be done by 
establishing relations between families of trajectory quadruplets and structures of periodic orbit pairs. For details see 
Appendix^] the results differ for the two universality classes. 

In the unitary case we find 

[O if to odd 1—1 if to odd . . 

dm = < , x m = i (40) 

II it to even , I (J it to even . 

The total contributions of all d- and x-families (per channel combination) now read 

D = — +V dm = 1 

TV 2 N m + 2 N 2 - 1 

ro+1 

oo ^ 

EsCfYi r ^ 

7V^+2 = ~N{N 2 - 1) ' (41 ^ 



m— 1 



The resulting conductance variance 



<(tr(ttt)) 2 ) _ ( tr(tt t)) 2 = i} (42) 



agrees with the random-matrix prediction l|32l) . 
In the orthogonal case we have 



om I I om _ i 

Xm = { _ ir ±__ (43) 



The contributions of d- and ir-families per channel combination now read 

TV 2 + 7v™+ 2 N(N + l)(N + 3) 



oo 

^JSfm+2 N(N +1)(N + 3) y ' 

and determine the variance in search as 

<(tr(«t)) 2) _ <tr(tt t))» = 2^(^ + 1)^ + 1) 

U V " 7 w ;/ 7V(iV+ l) 2 (iV + 3) ' V ' 

again in agreement with l|32|) . Thus, we have once more verified the universal behavior of individual chaotic cavities. 

IV. SHOT NOISE 

Our reasoning can be extended to a huge class of observables which are quartic in the transmission amplitudes and 
thus determined by d- and ^-quadruplets as well. For a first example, we consider shot noise: Due to the discreteness 
of the elementary charge, the current flowing through a mesoscopic cavity fluctuates in time as I(t) = I + SI(t) 
where / denotes the average current. These current fluctuations, the so-called shot noise, remain in place even at zero 
temperature. They are usually characterized through the power 



P = 4 / 6I(t )5I(t + t)dt (46) 
Jo 

where the overline indicates an average over the reference time to- 8 



This definition, as wel l as the treatment of three- lead correlation in the following section, follows the conventions of ITl l2fl . and differs 
by a factor 2 from . 
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Using our semiclassical techniques, we proceed to showing that for chaotic cavities, the energy-averaged power of 
shot noise takes a universal form. Again, our treatment applies to individual cavities and yields an expansion to all 
orders in the inverse number of channels. That expansion turns out convergent and summable to a simple expression 
which subsequently to our prediction was checked to agree with random-matrix theory by Savin and Sommers psj . 

Following Biittiker [l^ , we express the power of shot noise through the transition matrices £ 

(P) = (tr(«t) - tr(ttt tt t)) ; (47) 

here, P is averaged over the energy and measured in units 2e J^ depending on the voltage V. While the average 
conductance (tr(tft)) was already evaluated in Section^ the quartic term turns into a quadruple sum over trajectories 
similar to the conductance variance 

<(tr(ttt tt t))) = ( £ f aia2 t* l02 W: iC2 ) = ( E E AvA^A^-^-^A ; (48) 

- c 2 /3:ci^a 2 
(5: a\^C2 



the trajectories a, (3, 7, S must now connect the ingoing channels a\,c\ to the outgoing channels 02,02 as indicated 
in the summation prescription. 

As a consequence, the possible channel combinations for d- and ^-families of quadruplets are changed relative to 
the conductance variance. In the present case, d- quadruplets, with one partner trajectory coinciding at its beginning 
and end with a, and the other partner trajectory doing the same with 7, contribute only if either the ingoing or the 
outgoing channels coincide. If the ingoing channels coincide, the partner trajectory connecting the same points as a 
is of the type a\ = c\ — > 0,2 and may be taken as (3, whereas the trajectory connecting the same points as 7 has the 
form a\ = c\ — ► C2 and may be chosen as 5. If the outgoing channels coincide, similar arguments apply, with (3 and 8 
interchanged. Thus, ^-quadruplets contribute only for i\qiV| + N1N2 = NN1N2 channel combinations. In this sense, 
they take the role played by ^-quadruplets in case of the conductance variance. 

In turn, x- quadruplets now contribute for all channel combinations. Moreover, if both the ingoing and the outgoing 
channels coincide, either of the two partner trajectories may be chosen as /3 or S, meaning that the corresponding 
channel combinations have to be counted for a second time. Thus, ^-quadruplets now contribute for altogether 
N1N2 + N1N2 channel combinations, like ci-quadruplets in case of the conductance variance. 

We can simply interchange the multiplicity factors in our formula for ((tr(tt"f)) 2 ), Eq. gg), to get 

(tr(ttHt r )) = NN X N 2 D + (iV-fiVf + N 1 N 2 )X (49) 

and thus 

/pv _ { N N (I$-i) unitary case 

\ r l - \ N,(N,+l)N,(No + l) ., \ 0V ) 



N 1 (N 1 + 1)N 2 (N2 + 1 
N(N+l){N+3) 



orthogonal case 



Eq. (|50(l extends the known random-matrix result 

AT? AT. 



O ( ) unitary case 



N£Nj , N 1 N 2 (N 1 -N 2 ) 

"Tv 3 1 W 



O^j^j orthogonal case , 



to all orders in j?, for individual chaotic cavities. We can, moreover, give an intuitive interpretation for the terms 
in l|51|l . The diagonal contributions to (tr(££ T )) and {tr(ttHtf)) both read N ]^ 2 and therefore mutually cancel. The 

leading contribution, j^ 2 , arises from <i-quadruplets differing in a single 2-encounter (see Fig. 0J). In the unitary 
case, there are no terms of order 1, since all related families require time- reversal invariance. In the orthogonal case, 
Richter/Sieber pairs yield a contribution — N ^ 2 to (tr(fft)), from which we have to subtract two contributions to 
(tr (ttHtf ')) , the term — 2N ^ 2 accounting for (^-quadruplets differing in a single antiparallel 2-encounter (see Fig. 0Jd), 

and a term ^^f 2 arising from ^-quadruplets contributing to X2 = 4. The latter x-quadruplets may differ in two 
2-encounters, as in Figs. and [SJd, or in one 3-encounter, as in Fig. [SJ;. From the examples in Fig. [5] further 
families are obtained by interchanging a and 7, interchanging the two leads (for Figs. and[S|;), or interchanging 
the pairs (a, 7) and {(3,8) (for Fig. |SJd). Each of the Figs. |3Ji-c therefore represents altogether four families, whose 




FIG. 5: Families of ^-quadruplets with m — L — V = 2, contributing to the next-to-leading order of shot noise for time-reversal 
invariant systems. 




FIG. 6: A cavity with three leads. If a voltage V is applied between lead 1 and leads 2 and 3, one observes currents and 
/C 1 ^ 3 ) a s explained in the text, correlations between these currents are again determined by families of d- and ^-quadruplets 
of trajectories. 



contributions indeed sum up to = 4(— 1) 2 + 4(— 1) 2 + 4(— 1) = 4. Together with the contributions mentioned before, 
they combine to H51|l. 9 



V. CURRENT CORRELATIONS IN CAVITIES WITH THREE LEADS 



Another interesting experimental setting involves a chaotic cavity with three leads, respectively supporting iVi , N2 
and N3 channels; see Fig. [S] The second and the third lead are kept at the same potential, and a voltage is applied 
between these leads and the first one. Consequently, currents /( 1_>2 ), /t 1 ^ 3 ) flow from the first lead to the second 
and third one. We shall be interested in the fluctuations SI^ 1 ^ 3 ^ of these currents around the corresponding 

averages values, and study correlations between W^ 1- * 2 ) and W^ 1- * 3 ) [2^. l30| . 

This setting is similar to the famous Hanbury Brown-Twiss experiment [31| in quantum optics: there, light from 
some source (corresponding to the first lead) was detected by two photomultipliers (corresponding to the second and 
third lead) . Similar work on Fermions began somewhat later 29, 30] but the precise from of the correlation function 
is as yet unknown. 

Our semiclassical reasoning can easily be extended to fill this gap. The two currents depend on the matrices t^-^ 2 \ 
£(i->3) containing the transition amplitudes between channels of the first and the second and third lead; these matrices 
have the sizes N\ x N2 and N\ x JV3. As shown in |20L l30| . correlations between 5/( 1- * 2 ) and SI^ 1 ^^ are determined 
by the transition amplitudes as 

4 r 5I(^) (t )Sl(^) (t + t)dt = -(tr^^^^V 1 ^ 1 ^)) (52) 
Jo 



In ll.l . trajectory quadruplets where the encounter directly touches the lead are shown to become relevant when the mean dwell time 
is of the order of the Ehrenfest time. 
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(in units of a )• Using the semiclassical expression for the transition amplitudes, we are again led to a sum over 
quadruplets of trajectories 



(tr^^V^V 1 ^)) = ( E 



ai ,ci — 1.. . Ni 
a 2 = l...N 2 , 
c 3 = l...N 3 



-L( E E A a A^ 7 A|eH^+^)/A , (53) 



T 



\ai,ci,a 2 ,C3 a:ai^a 2 
/3: ci^a 2 
7: ci^c 3 
(5: ai — >C3 



with ai, ci, 02, C3 labelling channels of the first, second and third lead, as indicated by the subscript. The trajectories 
a, f3, 7, S must connect these channels as a(a\ — ► 02), /3(ci — > 02), 7(ci — > C3), <5(ai — > C3). 

The contribution of each family of trajectory quadruplets can be evaluated similarly to the conductance variance 
or shot noise. Since a particle can leave the cavity through any of the three leads, the escape rate depends on the 
overall number of channels iV = TVi + 7V2 + -/V3. Again, integration brings about factors i and —N for each link and 
each encounter. Only the numbers of channel combinations are changed, x- quadruplets as in Fig. [SJi contribute for all 
N1N2N3 possible choices of a%, ci, 02, C3. For any of these choices, partner trajectories connecting the initial point of 
7(01 — > C3) to the final point of a(ai — > 02), and the initial point of a to the final point of 7 are of the form ci — * 02 
and ai — ► C3 and can be chosen as f3 and S, respectively. In contrast, d- quadruplets as in Fig. Eh contribute only for 
the N1N2N3 combinations with coinciding ingoing channels a\ = c\. For these combinations the partner trajectory 
coinciding at its ends with a is of the type a\ = c\ —* 02 and can be taken as /3 whereas the partner trajectory 
coinciding at its ends with 7 has the form a% = c% — * C3 and can be chosen as 5. With these numbers of channel 
combinations, the current correlations in a 3-lead geometry are obtained as 



VI. ERICSON FLUCTUATIONS 



Another interesting quantum signature of chaos are so-called Ericson fluctuations, which have first been discovered 
experimentally in nuclear physics. In compound-nucleus reactions with strongly overlapping resonances, universal 
fluctuations in the correlation of two scattering cross sections at different energies have been observed. A first 
interpretation in terms of random-matrix theory was provided by Ericson |32| and further theoretical investigations 
have been reported in |3^|, [34[. 

Later on, the relation between classical chaotic scattering and Ericson fluctuations in single-particle quantum 
mechanics has been discussed j3^- Theoretical work shows that, e.g. the photoionization cross section of Rydberg 
atoms in external fields show universal correlations once the underlying classical dynamics is chaotic |3|J . 

Chaotic transport through a ballistic cavity displays Ericson fluctuations in the covariance of the conductance at 
two different energies, 



(C(E,e))^(^G(E)G^E+^j^-(G) 2 , (55) 

with G(E) = tr(tv). Here, the difference between the two energies was made dimensionless by referral to the energy 
scale 2^ proportional to the number of channels and to the mean level spacing. Similarly as for the conductance 

variance, the semiclassical approximation for t(E),t(E') where E' = E + leads to a quadruple sum over 
trajectories, 

(tvtt^E) tvtt\E')) = ( E t aia AE)tl iaa (E)t ClC2 (E')t: iCa (E')) (56) 



«2 > c 2 



= J_/ V E A a A*pA 1 Ale l[Sa{E) - Sl){E)+s -< {E ' ) - Ss( ' E ' ))/n 

H \ 01,01 a f3; ai ^a 2 I 
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the only difference to l|33[) being that the trajectories 7 and 5 have to be taken at energy E' = E + ^L. Using 
= T 7 , = Ts , and Th — lufip, the phase factor can be cast into the form 

e i(S a (E)-S fi (E)+S^(E+^)^S 6 (E+^))/h _ e i(S a (E)-Sf,(E)+S^(E)-S s (E))/H x e(T 7 -T s ) ^ ^ 

i.e., the quadruple sum in (|56(l differs from l|33|l by an additional factor depending on the difference between the dwell 
times of 7 and S. 

The latter difference may be written as a sum over links (with durations £j) and encounters (with durations ienc)i 

L+2 V 

T 1 -T s = J2viti + Y, T l°%n C - (58) 

1=1 (T=l 

Here, the integer numbers r\i and rj a characterize the individual links and encounters. (Note the distinction between 
links and encounters by Latin and Greek subscripts). Each link occurs twice in the quadruplet, once in one of the 
original trajectories a, 7 and then in one of the partner trajectories /3, S. The number r\i = 0, ±1 gives the difference 
between the numbers of times the i-th link is traversed by the trajectories 7 and S. We thus have rji = 1 if the i-th 
link is traversed by 7 and not by 5, i = if it is traversed either by both or none of the two trajectories, and i = — 1 
if it is traversed only by 5. Similarly, r\ a gives the difference between the numbers of traversals of the cr-th encounter 
by 7 and 5. For an ^-encounter, r\ a may range between —I and I. 

When evaluating the contribution of each family of quadruplets, Eq. we simply have to add a phase factor 

e t h " 1 for each link, and a phase factor e T « onc for each encounter. The link and encounter integrals are thus 
replaced by 

at t e e ~ jv(1 - iry i£ ) ' 1 ' 



and 



/ 



SV-Hz nc (s,u) I T~ V ' 



As a consequence, our diagrammatic rules are modified to yield a factor jfjrh^r^. for each link and a factor —N(l — 
irja-e) for each encounter. 

We must, however, be aware that the numbers rji, r\ a depend on which of the two partner trajectories is labelled as (3 
and which is labelled as <5. Each family of d- or x-quadruplets hence comes with two different sets of numbers {rji, r\a\, 
depending on the combinations of channels considered. For each ci-family we have to keep into account NfN^ channel 
combinations with <5 coinciding at its ends with 7; all these "^7-type" combinations give rise to the same {?7i,?7o-} and 
to the same link and encounter factors. In addition, we must consider N1N2 combinations of the "fo-type" with S 
coinciding at its ends with a, and a different set of {77^,77^}. 

For each x-family we would, in principle, have to distinguish between N\N% combinations with coinciding ingoing 
channels, and 6 coinciding at its beginning with a and at its end with 7, and N^N2 combinations with coinciding 
outgoing channels, and S coinciding at its beginning with 7 and at its end with a. Such caution is, however, unnecessary 
for reasons of symmetry. Each x-family is accompanied by another one which is topologically mirror-symmetrical, 
with left and right in Fig. ^interchanged. In this family, initial points turn into final ones, and vice versa, implying 
that (3 and 5 are interchanged. Since both families are taken into account simultaneously, "mistakes" like always 
choosing S to connect the initial point of a to the final point of 7, are automatically compensated. 

We can thus write the conductance covariance as 

C(e) = N?Nl £ ^_ + Nl N 2 £ ^ + N lN2 N £ ^ - (G) 2 (61) 

m=0 m—0 m—1 

Here the coefficients x m (e),drn' 1 \e),drn a \e) are the summary contributions of the ^-quadruplets and the two men- 
tioned groups of (i-quadruplets, with m = L(v) — V(v) — 2 (and thus m — for the diagonal quadruplets) and the 
denominator N~ m ~ 2 dropped. The squared averaged conductance is e-independent and is determined by Eq. (J2J. 
The leading contribution to the conductance covariance corresponds to dropping in (|(jl|) all coefficients but 

x%, d^ Sl \ dtj Sa \ d£ Sl \ d^ 5 ^ ■ For the conductance variance, we had seen that the contributions of (^-quadruplets 
that fall into pairs (a, ft) and (7, S) contributing to conductance cancel with the squared average conductance. The 
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+i 

FIG. 7: Families of trajectory quadruplets contributing to the covariance of conductance (coinciding with Fig. EJi-hj). The 
trajectories 7 and 5 are highlighted through dashing and dotting, assuming that 5 connects the same points as 7. The picture 
also indicates all non- vanishing numbers rji and r] a (the latter in bold font). 



same remains valid here, since for these pairs 7 and S traverse the same links and encounters, and all rji and r\ a vanish. 
Again, the contributions of diagonal quadruplets a = S, ft = 7 and x-quadruplets as in Fig. (see also Fig. 0*) 
mutually compensate; compared to the variance both kinds of quadruplets receive the same additional factors due to 
one link with rji = 1 (the trajectory /3 = 7 and the lower left link in Fig. [7^) and one link with rji = — 1 (the trajectory 
a = S and the upper left link in Fig. 0j) . 

Like for the conductance variance, the leading contribution to the covariance thus originates from the d-families in 
Fig. 0J>h which contain encounters between a and 7 and thus do not fall into pairs relevant for conductance. These 
families are redrawn in Fig. [7^,-f, together with all non- vanishing numbers rji and rj a . The trajectories 7 and 6 are 
highlighted through dashing and dotting, assuming that S connects the same points as 7; channel combinations with 
S connecting the same points as a only contribute to higher orders in -k. The family of Fig.0L involves a link with 

r\i = — 1 (the upper central one) and a link with rji = 1 (immediately below), and thus yields ^e^i^^i_ ie -j — ' 
The same holds for the family in Fig. 0i, which requires time-reversal invariance. In contrast, the contributions of 
Fig. 03,c,e, and f remain independent of e, since additional factors from links with rji = — 1 and encounters with 
r/a- = —1 mutually compensate; the same applies for the families represented by "x2" in Fig. [7| with (a, (3) and 
(7, S) interchanged and the signs of r\i and r\ a flipped. As for the variance of conductance, the contributions of 
Fig. Et>,c,e,f thus mutually cancel, both in the orthogonal case and in the unitary case (where only Figs. Et>,c may 
exist). Altogether, we now obtain 



/ eN 

g(E)g(e + — 



jV<t(l+e 2 ) + U \n) 

2N ? N Z_ 1 n (±.\ 



unitary case 
orthogonal case . 



(62) 



The Lorentzian form of (|62|l confirms the random- matrix predictions of |3^| . 

Higher orders in , not known from random- matrix theory, can be accessed by straightforward computer-assisted 
counting of families of quadruples differing in a larger number of encounters, or in encounters with more stretches. To 
do so, we generated permutations which describe possible structures of orbit pairs (see Appendix iBl and Q). We then 
"cut" through these pairs as described in Appendices lAl and IBl to obtain quadruplets of trajectories and determined 
the corresponding rji, rj a . The final result can be written as 



n?n* r 1 

W 4 (l+e 2 ) ^ 
4N1N2 



_ 1 l+3£^+21e 4 +5e"+2e 8 I 
) 1 N^(l+^f / 

2JV1JV2 _ 2jV 1 2 JV 2 2 (5+12e 2 +3e 4 
Af 5 (l+e 2 ) a 



JV 4 (1 

+0(w 



Af 3 (l+e 2 ) 

2AfjAf 2 (l8+78e 2 + 177e 4 +48e 6 + llc 8 ) 



unitary case; 



orthogonal case . 



(63) 



20 



In the unitary case the a;- type contribution cancels in all orders with the d ^ Sa ^ -contribution; for that reason the overall 
result is proportional to NfN%. 

VII. QUANTUM TRANSPORT IN THE PRESENCE OF A WEAK MAGNETIC FIELD 

A. Changed diagrammatic rules 

Our methods can also be applied to the case of a weak magnetic field, with a magnetic action of the order of fi. 
The necessary modifications were introduced in [^3] for the spectral form factor; see also [2^. As in [23], we will 
obtain results interpolating between the orthogonal case (without a magnetic field) and the unitary case, where the 
magnetic field is strong enough to fully break time-reversal invariance. We shall assume that the field is too weak 
to influence the classical motion, meaning that we have to deal with the same families of trajectory pairs as in the 
orthogonal case. However, the action of each trajectory is increased by an amount proportional to the integral of the 
vector potential A along that trajectory, e.g., by 

Q = /-A(q)-dq, (64) 

for the trajectory a. When we evaluate the average conductance, the action difference inside each pair of trajectories 
a and (3 is thus increased by Q — Qp. This additional term may be neglected for pairs of trajectories where all 
encounters are parallel. For these pairs, all links and stretches of (3 are close in phase space to links and stretches of 
a, and therefore receive almost the same magnetic action. 

The situation is different for pairs where a and (3 traverse links or stretches with opposite sense of motion. Since 
the magnetic action changes sign under time reversal, such orbit pairs have significant magnetic action differences 
@q — 0/3- These differences can be split into contributions from the individual links and encounters. Let us first 
consider links. If (3 contains the time-reversed of the i-th link of a, it must obtain the negative of the corresponding 
magnetic action 0^. The difference Q — 0/3 then receives a contribution 20^. Therefore we may write the contribution 
of each link as 2/i^0i with \ii = 1 if the link changes direction on (3 and [ii = otherwise. 

Consider now the contribution of encounters. We assume that in the original trajectory a the encounter a had v a 
stretches traversed in some direction (arbitrarily chosen as "positive" ) meaning that the remaining l a — v a stretches 
were traversed in the opposite, "negative" direction; in the trajectory partner (3 these numbers will generally change 
to v' a ,la — v' a correspondingly. Denoting the magnetic action accumulated on a single stretch traversed in a positive 
direction by CT we see that the encounter a yields 2/v0<r to the magnetic action difference, with \i„ = v' a — v a . The 
overall magnetic contribution to the action difference now reads 

L+l V 

Q - ©/s = 2 ^ e * + 12 2 ^ e ° ( 65 ) 

1=1 (7=1 

and yields a phase factor 

L + l V 

J] e l2 ^' h Y[ e i2 ^ e ^ h , (66) 

1 = 1 (T = l 

where we again distinguish between links and encounters only through Latin vs. Greek subscripts. 

To handle this additional phase factor, we show that for fully chaotic (in particular, ergodic and mixing) dynamics, 
the magnetic action may effectively be seen as a random variable [37|. For fully chaotic systems, any point on any 
trajectory can be located everywhere on the energy shell, with a uniform probability given by the Liouville measure. 
Moreover, phase-space points following each other after times larger than a certain classical "equilibration" time t c \ 
can be seen as uncorrelated. We will therefore split each link or encounter stretch into pieces of duration t c \. These 
pieces have different magnetic actions. Let us consider the probability density for these actions. Since positive and 
negative contributions to the magnetic action are equally likely, the expectation value for the action of an orbit piece 
must be equal to zero. The width W (i.e., the square root of the variance W 2 ) must be proportional to the vector 
potential and therefore to the magnetic field B. Since the magnetic actions of the individual pieces are uncorrelated, 
the central limit theorem then implies that the magnetic actions of links with K = 3> 1 pieces obey a Gaussian 

probability distribution with the width \f~KW , i.e., 
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The phase factor arising from a link averages to J dQiP{Qi)e l2 ^ iBi ' h = e ^ ibti : depending on the system-specific 

parameter b = ^faj— = jp|- oc -p- and on /Zj = /zf G {0, 1}. Similarly, the phase factor associated with the er-th 

encounter averages to e~^ 6t ~"= 37]. Links and stretches traversed in opposite directions by a and (3 thus lead to 
exponential suppression factors in the contributions of trajectory pairs. 

These factors have to be taken into account when evaluating the average conductance, starting from 122fl . The link 
integrals are changed into 

> N rp 

= I , (68) 
with £ = -^-b oc tU, whereas for each encounter we find an integral 

\7 ^-i^ nc ( S , M ) / T~ v 7 

Since the Tjf's again mutually cancel, our diagrammatic rules are changed to give a factor jv(i+/j for each link 
and a factor —N(l + fi^C) fo r each encounter; the arising product has to be multiplied with the number of channel 
combinations, i.e., N1N2 for the average conductance. 

The same rules carry over to the conductance variance, shot noise, and correlations in a three-lead geometry. In 
these cases, /x, is equal to 1 if the i-th link of the pair (a, 7) is reverted in (/3, S), and counts the stretches of the 
CT-th encounter of (a, 7) which are reverted in (/3, 6); the sign of \i a is fixed as above. 



B. Mean conductance 

For the average conductance, the diagonal contribution, Nl j^ 2 , remains unaffected by the magnetic field. The 
contribution of Richter/Sieber pairs, — N ^ 2 , obtains an additional factor j^, since one of the three links of a in 
Fig. nor|2t is traversed by f3 in opposite sense. The next order originates from trajectory pairs as in Fig. [2t>-j where 
arrows indicate the direction of motion inside the encounters and highlight those links which are traversed by a and (3 
with opposite sense of motion. The contributions of the families in Fig. [2b, c, d, g, i, j remain unchanged: In Fig. [2t>, 
g no links or encounter stretches are reverted; for Fig. Et, d, i, j the number of links with yn = 1 and encounters with 
//J. = 1 coincide, meaning that the ^-dependent factors mutually compensate. The six above families cancel mutually 
due to the negative sign for Fig. i, j- The contributions of Fig. [5^, f, h obtain a factor mjgp from two reverted 

links; due to the negative sign of Fig. |2h, they sum up to ^^^2 ■ We thus find 

M«'» = ^{i-^ + w ^} + °(>)- (*» 

Counting further families of trajectory pairs with the help of a computer program one is able to proceed to rather 
high orders in ■h . We then find 

1 1 1 l + 2£ + 13£ 2 + 4£ 3 +£ 4 1 l + 2£ + 49£ 2 + 4£ 3 +£ 4 

AT 2 " (1 + 2 ~ AT 5 (1 + 5 + AT* (1 + C) 6 

(71) 

As expected, (I7UII and l|71l) interpolate between the results for the orthogonal case, reached for B — * and thus 
£ — * 0, and the unitary case, formally reached for B — * 00 and thus £ — ► 00. The convergence to the unitary result 
is non-trivial: The contributions of the families in Fig. |2t,d,i,j are not affected by a magnetic field, because all in- 
dependent factors cancel. These contributions thus survive in the limit £ — > 00 (i.e., when the magnetic action becomes 
much larger than h, but the trajectory deformations due to Lorentz force can still be disregarded), but vanish in the 
unitary case (i.e., when the magnetic field is strong enough to considerably deform the trajectories). The agreement 
between the limit (-*oo and the unitary result implies that the contributions of all such families must sum to zero, 
for all orders in -A.. Order by order, l|71|) coincides with the results of |3^|, where the individual coefficients were given 
as (rather involved) random-matrix integrals. 



tr (ttf) ) = — —I 1 - — - + 

K y " N y AT (1+0 
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C. Conductance variance, shot noise and three-lead correlations 



For observables determined by families of trajectory quadruplets, it is convenient to first evaluate the overall 
contributions of d- and rr-families per channel combination. These contributions, denoted by D and X, now depend 
on the parameter £. The contribution of d-families reads 



D 



1 

7p ~ 
1 

+ AT5 



AT3(l + £) 
10 



1 

W 

4 



1 



(i + e) 2 . 

2(1 + 4^ 



2(1 + 9Q 

i + £ (i + 3 (i + 5 "~ (i + 4 



(72) 



generalizing our previous results 141|) and l|44|l for the unitary and orthogonal cases. The leading term, originating 
from diagonal quadruplets, remains unaffected by the magnetic field. The second term is due to quadruplets as in 
Fig- Wp- Since in these quadruplets, one link of (a, 7) is time-reversed in (/3, 5), the corresponding contribution is 
proportional to j^. It is easy to check that the third term correctly accounts for <i-quadruplets differing in two 2- 
encounters, or in one 3-encounter; compare Subsection IIII Bl and Fig.0] The higher-order terms were again generated 
by a computer program. 

In the overall contribution of x-families, 



X = 



1 



1 

AT5 



-1 - 



10 



2(l + 4fl 

(i + 2 (i + O 4 



(73) 



the term —-rm accounts for x-quadruplets differing in a parallel 2-encounter, Fig. 0]. These quadruplets are not 
affected by the magnetic field. All families responsible for the second term, Fig. EJi-c, display a Lorentzian field 
dependence j^: While Figs. and|SJ: contain one time-reversed link and only encounters with fx a = 0, Fig. 03 

involves two time-reversed links and one encounter with [i 2 = 1 . The remaining terms were again found with the help 
of a computer. 

With these values of D and X, we obtain, writing out only terms up to 0(N^ 1 ), 



the conductance 



variance 



((tr(ttt)) 2 )_(t r(tt t)) 2 = ^l 



1 



(i + C) 5 



2N?N$ (5 + 8£ + 4£ 2 ) 2Ni 7V 2 



iv 5 (i + 3 



iV3(l+0 



o 



1 

N 2 
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• the power of shot noise 

NfN$ 7Vi7V 2 (7Vi - N 2 ) 2 N?N$(13 + 32£ + 160 + 40 + £ 4 ) 



(tr(tt t )-tr(ti t tt t )) - 



N 3 



7V 4 (l + 



• and current correlations for a cavity with three leads 



37ViiV2 
> 3 (l + 2 



/ tr(t (i^2) t (i^2)t t (i^3) t (i^ 3 )tjv = NiN 2 N 3 (N 2 + N 3 ) , 2N.NMN, - N 2 - N 3 ) 



N 3 7V 4 (l + 

TViAWs [(N 2 + JV 3 )(1 + 2 (5 + 2g + Q) - 2JVi(4 + lOg + 3Q)] /J_ 

7V 5 (l + 4 + \N 2 



O ^ 



1 

N 2 



(75) 



(76) 



At least the higher orders in j? are new results. In particular, for the power of shot noise, we do not only obtain 
the previously known cancellation of the second term at N\ — N 2 — y , but also a new field dependence due to the 
third term, 



8g + 4Q + 4Q + Q 
16(1 + 4 



O 



(n 2 ) ; 



(77) 



D. Ericson fluctuations 



When studying Ericson fluctuations in a weak magnetic field, we have to deal with two parameters (apart from the 
channel numbers): the scaled energy difference e and the parameter £ proportional to the squared magnetic field. Our 
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Contribution of each link 


simplest case 


i 

N 


with energy diff. cx e and squ. magn. field oc £ 


1 


Contribution of each encounter 


simplest case 


—TV 


with energy diff. oc e and squ. magn. field oc £ 




Total contribution 

per channel combination 


trajectory pairs 


unitary case 


i 

JV 


orthogonal case 


i 

N+l 


d-quadruplets 


unitary case 


1 


orthogonal case 


N+2 
N(N+l)(N+3) 


^-quadruplets 


unitary case 


1 

N(N 2 -1) 


orthogonal case 


AT(JV+l)(JV+3) 


Number of channel combinations 


trajectory pairs 


conductance 


NiN 2 


d-quadruplets 


variance of conductance 


N?N$ + 7Vi7V 2 


shot noise 


N!N 2 N 


3-lead correlations 


NxNiNs 


^-quadruplets 


variance of conductance 


N!N 2 N 


shot noise 


N?N% + NiN 2 


3-lead correlations 


N?N 2 N 3 



TABLE I: Diagrammatic rules determining chaotic quantum transport. The table shows link and encounter contributions for 
all observables discussed in the present paper. For the simplest cases (conductance, conductance variance, shot noise, 3-lead 
correlations), we have also listed the summed-up contributions of trajectory pairs and d- and x-quadruplets and the numbers 
of channel combinations. 



diagrammatic rules are then changed in a straightforward way. Each link yields a factor j^j^^z^ ^ , whereas each 
encounter gives — N(l + fi a S, — i^o-e), to be multiplied with the number of channel combinations. 

As in the orthogonal and unitary cases, the leading contribution can be attributed to the quadruplets in Figs. [7^ 
and0i; all other contributions of the same or lower order, including the remaining families in Fig.0 mutually cancel. 

__ N 2 N 2 

Quadruplets as in Fig. [Tp, do not feel the magnetic field, and thus yield j^r^ri as shown in Section IVII here we 
dropped lower-order corrections due to the case of coinciding channels. For the family of quadruplets depicted in 
Fig. EJi, the two links with r/i — ±1 connecting the two encounters are reversed inside (/3, S) and thus have /ii = 1. We 
obtain factors Af ^ 1+ 1 ^_ ie ^ and N ^^ +ie ^ from these two links, from each of the four remaining link, and — N from 

N 2 N 2 

the encounter. Multiplication with the number of channel combination yields a contribution N 4.^ 1 ^^y2 +e i^ ) ■ Ericson 
fluctuations in a weak magnetic held are therefore determined as 

( G(E)G ( E + - ) ) _ <G) . .HM + _L_ } + ( i ) . 



VIII. CONCLUSIONS 



A semiclassical approach to transport through chaotic cavities is established. We calculate mean and variance of 
the conductance, the power of shot noise, current fluctuations in cavities with three leads, and the covariance of the 
conductance at two different energies. These observables are dealt with for systems with and without a magnetic field 
breaking time-reversal invariance, as well as in the crossover between these scenarios caused by a weak magnetic field 
leading to a magnetic action of the order of H. In contrast to random-matrix theory, our results apply to individual 
chaotic cavities, and do not require any averaging over ensembles of systems. Moreover, we go to all orders in the 
inverse number of channels. 

Transport properties are expressed as sums over pairs or quadruplets of classical trajectories. These sums draw 
systematic contributions from pairs and quadruplets whose members differ by their connections in close encounters, 
and almost coincide in the intervening links. The contributions arising from the topologically different families of 
quadruplets or pairs are evaluated using simple and general diagrammatic rules, summarized in Tab. [I] (These rules 
remain in place even for observables involving higher powers of the transition matrix, as shown in Appendix lOl . 
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Our work shows that, under a set of conditions, individual chaotic systems demonstrate transport properties devoid 
of any system-specific features and coinciding with the RMT predictions. An obvious next stage would be investigation 
of the system-specific deviations from RMT observed when these conditions are not met. Previous work PIEIGI has 
already motivated an extension to the regime where the average dwell time To is of the order of the Ehrenfest time Te, 
i.e. the duration of the relevant encounters. Here, the semiclassical approach helped to settle questions controversial 
in the random-matrix literature . As shown in |l2l , the leading contributions to the average conductance and 
the power of shot noise become proportional to powers of e~ TE / Tr> , ultimately arising from the exponential decay of 
the survival probability. On the other hand, the conductance variance turned out to be independent of Te Il2l . 

The door is open for a semiclassical treatment of many more transport phenomena, such as quantum decay |40j , weak 
antilocalization |4l|, parametric correlations the full counting statistics of two-port cavities, and cavities with 

more leads. Extensions to the symplectic symmetry class, along the lines of 0,^3] j an d to the seven new symmetry 
classes [44| (relevant e.g. for normal- metal/superconductor heterostructures or quantum chromodynamics) should be 
within reach. A generalization to quasi one-dimensional wires would finally lead to a semiclassical understanding of 
dynamical localization. 

We are indebted to Dmitry Savin and Hans-Jiirgen Sommers (who have reproduced our prediction (|5L)|) in random- 
matrix theory 28]); to Piet Brouwer, Phillippe Jacquod, Saar Rahav, and Robert Whitney for friendly correspondence; 
to Taro Nagao, Alexander Altland, Ben Simons, Peter Silvestrov, and Martin Zirnbauer for useful discussions; to 
Austen Lamacraft for pointing us to 2.7j; and to the Sonderforschungsbereich SFB/TR12 of the Deutsche Forschungs- 
gemeinschaft and to the EPSRC for financial support. 



APPENDIX A: TRAJECTORY QUADRUPLETS VS ORBIT PAIRS 

In this Appendix we will establish a combinatorial method for counting families of trajectory quadruplets appearing 
in the theory of conductance variance and shot noise. We will see that trajectory quadruplets can be glued together 
to form orbit pairs, and orbit pairs can be cut into quadruplets of trajectories. In contrast to the case of trajectory 
pairs, see Fig. |3 we shall now need two cuts. 

Our approach will be purely topological; e.g., an orbit pair (A,B) is regarded just as a pair of directed closed 
lines with links coinciding in A and B but differently connected in the encounters. Similarly, within each quadruplet 
(a,/3, 7,5) we can assume that the links of a, 7 exactly coincide with those of (3,8. Mostly, we can even think of the 
quadruplets as black boxes with two left ports oi, c\ and and two right ports a%, ci- Regardless of the actual number 
of encounters inside, an x-quadruplet can then be treated like a "dressed" 2-encounter: connections a\ — a2, c\ — C2 in 
one of the trajectory pairs are replaced by a\ — C2, c\ — ai in the partner pair, exactly as if a single 2-encounter existed 
between the trajectories a, 7 of the quadruplet. On the other hand, no change in the connections occurs between the 
ports in a d-quadruplet, hence it is topologically equivalent to a pair of dressed links. 

We shall consider both the unitary and the orthogonal case. In each case, we will use two slightly different methods 
to relate trajectory quadruplets and orbit pairs. This will allow us to express the quantities x m and d m defined in 
(|37|l through the auxiliary sums 

M(v)=m 

A m = £ (-l) v ^(L(v) + l)Af(v) 

V 

M(v)=m 

Bra = £ (~l) y( ^(«,2), (Al) 

v 

where Af(v) and Af(v, 2) are numbers of structures of orbit pairs (see Subsection lll D|) and we have M(v) = L(v) — V(v); 
these auxiliary sums will be determined recursively in Subsection IA 31 below. 



1. Unitary case 

To illustrate method 7, let us consider a d-quadruplet (a,/3, 7, 5), as on the left-hand side of Fig. [3] and merge a 
and 7 into one "orbit" A. We connect the final point of a to the initial point of 7 and the final point of 7 to the 
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FIG. 8: Left-hand side: Schematic picture of a d-quadruplet of trajectories a, 7 (full lines), /3, 8 (dashed lines). Right-hand 
side: Outside the "bubble", we added connection lines joining (0,7) and (/?, 5) into periodic orbits A and B with the initial 
link indicated by "1.". The additional lines of A and B coincide. 
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FIG. 9: Left-hand side: Schematic picture of an ^-quadruplet of trajectories 0,7 (full lines), /3, S (dashed lines). Right-hand 
side: Outside the "bubble", we added connection lines joining (0,7) and ((3,6) into periodic orbits A and B. The additional 
lines of A and B differ from each other, and can be viewed as an additional 2-encounter. 

initial point of a, as shown on the right-hand side. Likewise, /3 and S can be glued together to an "orbit" 23. 10 The 
connection lines added are the same for (a, 7) and for (/3, 5): one connection line joins the coinciding final links of a 
and /3 with the coinciding initial links of 7 and 5, whereas the second one joins the final links of 7 and 5 with the 
initial links of a and j3. The orbits A and B differ in the same encounters as (a, 7) and (/3, 5). To fix one structure for 
the orbit pair A, B, we have to single out one link as the "first" and choose as such the link of A created by merging 
the final link of 7 with the initial link of a (indicated by "1." in Fig. |SJ. 

We can revert the above procedure, to obtain families of d-quadruplets from structures of orbit pairs. We first have 
to cut both orbits inside the "initial" link. This leads to a trajectory pair with L(v) + 1 rather than L(v) links. We 
then have L(v) + 1 choices for placing a second cut in any of these links. In each case, we end up with a trajectory 
quadruplet. Within this quadruplet, the trajectories following the first cut through A and B are labelled by a and 
(3; the remaining ones are called 7 and S. In this way, each of the Af(v) structures of orbit pairs related to a given 
v gives rise to L(v) + 1 families of ^-quadruplets with the same v. The quantities Md{v) and d m characterizing the 
d-familics in (|37|) thus become accessible as 

Ki(v) = (L(v) + l)M(v) , (A2) 

L(v)-V(v)=m 

d m = Y, (-l) V(vl (m + lW(v) = A m . (A3) 

V 

Let us discuss a few examples. According to (IA3|) the coefficient d\ is determined by orbit pairs with m = 1 
(i.e., only one 2-encounter). Since in the unitary case there are no such orbit pairs, we have d\ = 0. The following 
coefficient di is determined by orbit pairs with m = 2. In the unitary case there are only two such structures, ppi and 
pc (Fig. [5Jd and g). All quadruplets responsible for the coefficient c?2 can be obtained by making two cuts through 
these orbit pairs, one through the initial link which may be chosen arbitrarily. The second cut can go through any 
link; in particular, there are two possibilities for the second cut in the initial link, before and after the first cut. That 
means L + 1=5 possible positions of the second cut for ppi. These lead to the quadruplets as in Fig. 2t and d, the 
reflected version of Fig. and quadruplets where either a or 7 contain two 2-encounters and the other trajectory 



As mentioned, f} and 6 may be interchanged if the ingoing and outgoing channels coincide. This has no impact on the present 
considerations. The naming of partner trajectories as /3 and 6 in all figures will be arbitrary. 
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FIG. 10: Sketches of orbit pairs (.4, B). Inside the bubbles: quadruplets of trajectories obtained by cutting (A, B) in the initial 
link (indicated by "1."), and in one further link. In (a) both links are traversed by A and B with the same sense of motion, 
and the trajectory quadruplet is of type d. In (b), the second cut is placed in a link traversed with opposite senses of motion, 
and the resulting quadruplet is of type x. 

contains none. For pc there are four possible positions for the second cut, corresponding to Fig. 0^, its reflected 
version, and quadruplets where either a or 7 contain the full 3-encounter. All quadruplet families related to a given 
structure make the same contributions (— \) v to the coefficient d 2 , i.e., 1 for those obtained from ppi and -1 for those 
obtained from pc; we again see that cfe = 5 — 4 = 1. 

To explain method II, let us now consider x- quadruplets as on the left-hand side in Fig. EI O n the right-hand side, 
a and 7 are again merged into a periodic orbit A, and (3 and 5 are once more merged into B, by connection lines 
leading from the end of one trajectory to the beginning of the other one. In contrast to the first scenario, the pair 
(A, B) has one further 2-encounter between these lines, with different connections for the two partner orbits. To fix 
one structure for the latter orbit pair, we take the initial link of a as the "first" link of the orbit pair. This link is 
preceded by a "final" stretch, which must belong to the added 2-encounter. We must therefore reckon with orbit pairs 
associated to the vector z/[~* 2 l and whose final stretches belong to a 2-encounter. In the notation of Subsection III Dl 
the number of structures of such orbit pairs is given by Af(v^ 2 ^ , 2). 

Each of these structures can be turned back into one family of x-quadruplets, by cutting out the added 2-encounter. 
Consequently, there is a one-to-one relation between x-families and the structures of orbit pairs considered. The 
number of ^-families related to v and the coefficients x m defined in (|37[1 are given by 

M x (v) = 2), (A4) 

L(v)-V{v)=m 

(-l) v( ^(^ a] ,2). (A5) 

V 

Rather than over v, we may sum over v' = t/^ 2 ! with L(v') — V(v') — (L(v) + 2) — (V(v) + 1) = m + 1. While the latter 
sum should be restricted to v' with v' 2 > 0, that restriction may be ignored since v' with v' 2 = have Af(v',2) = 0. 
Using (—l) v ( v > = — (— l) v ( v ) and dropping the primes we can express the coefficient x m as 

L(v)-V(v)=m+1 

x m = - {-l) V{0) M{v,2) = -B m+1 . (A6) 

v 

For instance, the coefficient x\ is determined by orbit pairs with L(v) — V(v) = 2. Only one such structure (ppi) 
contains 2-encounters and yields a contribution —1, whereas pc involves only one 3-encounter. We thus find x\ = —1, 
as already seen previously. 

Taken together, Eqs. (|A3I) . I|A6|) . and (|2*FJ)l indeed relate the conductance variance to structures of orbit pairs. 

2. Orthogonal case 

For time-reversal invariant systems, we must now consider pairs of orbits A, B differing in encounters whose stretches 
are either close or almost mutually time-reversed. The sense of traversal of an orbit now being arbitrary we fix the 
direction of B such that B traverses the "initial" link of A in the same direction. 

We start with method I, i.e., we cut an orbit pair (A, B) inside links, first inside the initial link and afterwards in 
an arbitrary link. We have to distinguish two cases, respectively leading to d- and x-quadruplets. First assume that 
the second cut is placed in a link traversed by A and B with the same sense of motion; see Fig. ISlor llUb . As in the 
unitary case we then obtain a d-quadruplet of trajectories; this quadruplet is highlighted by a grey "bubble" in Fig. [SI 
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FIG. 11: Sketches of orbit pairs (A,B) with one 2-encounter singled out. This 2-encounter contains the "final" encounter 
stretch, and is either (a) parallel in both A and B, (b) parallel in A and antiparallel in B, (c) antiparallel in A and parallel in 
B, or (d) antiparallel in both A and B. 



and HOfa . without depicting the encounters. Inside this quadruplet, a is defined as the trajectory following the cut 
through the initial link (marked by "1." in Fig. 1811 1|) . 

Now assume that the second cut is placed in a link traversed by A and B with opposite senses of motion. The 
resulting quadruplet, shown in the grey bubble in Fig. 110b . resembles an ^-quadruplet, apart from the directions of 
motion. To obtain a true x-quadruplet, one has to revert the directions of motion of two trajectories in Fig. llUb . such 
that all trajectories point in the same direction as a (i.e., the trajectory following the cut inside the initial link). 

We hence obtain the following relation between orbit pairs and trajectory quadruplets: By cutting inside links, each 
of the Af(v) structures of orbit pairs related to v can be turned into a family of trajectory quadruplets in L(v) + 1 
possible ways. Through such cuts, all Afd{v) (i-families and all Af x (v) x-families are obtained exactly once, since each 
d- or x- family could be inserted inside the bubble in Fig. 110b orllOb. respectively. We thus have 

(L(v) + l)Af(v) = Af d (v) +Af x (v) (A7) 

and summation as in the unitary case leads to 

= d m + x m (A8) 

with A m defined in IjAlfl and d m , x m defined in (|37|) . 

We now turn to method II (cutting inside 2- encounters). We consider orbit pairs A, B containing one more 2- 
encounter compared with the quadruplets in question, assuming that the final encounter stretch of the orbits belongs 
to the "added" 2-encounter; the number of these structures is J\f(v^ 2 \2). Cutting out the added 2-encounter will 
create all possible quadruplets associated to v, some of them, as we shall see, in several copies. 

The added 2-encounter can be either parallel or antiparallel in the orbit A as well as in its partner B. This leads to 
four different possibilities, depicted by arrows on white background in Fig. Illfa -d. First suppose that the encounter 
in question is antiparallel in both A and B, as in Fig. II 1H . This is possible only if the ports of this encounter are 
connected in the same way in A and B, up to time reversal; one can easily check that all other connections would lead 
to either A or B decomposing into several disjoint orbits. The connections outside the encounter, as depicted in the 
bubble in Fig. Illt l. thus resemble d-quadruplets and can be turned into true ci-quadruplets if we revert the sense of 
motion on some trajectories. 11 Any rf-family could be substituted for the bubble in Fig. I11H . Thus, cutting through 
2-encounters of the kind in Fig. Illt l produces all possible families of d-quadruplets. 



11 Rather than reverting directions of motion, we could also identify the initial points of all trajectories inside the bubble with ingoing 
leads, and the final points with outgoing leads, loosing the identification of the two sides of our bubble with the two openings of the 
cavity. This would entail a different mapping between orbit pairs and trajectory quadruplets, but not affect the following results. 
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In the three other cases, Fig. Illh-c. the remaining connections have to be of type x, up to the sense of motion on 
some trajectories, since connections of type d would lead to decomposing orbits. To better understand these cases, it 
is helpful to view the corresponding x-quadruplets as a "dressed" 2-encounters. Then, the orbit pairs of Fig. ITTk-c 
are topologically equivalent to the simple diagrams in Fig. [21 The orbit pair in Fig. Illfa is of the type ppi whereas 
the pairs in Fig. II lb and c are both of the type api. (In Fig. II lb the initially parallel encounter is identified with the 
added encounter, and the initially antiparallel encounter is identified with the ^-quadruplet, whereas the situation is 
opposite in Fig. 111b. ) Any family of type x could be substituted for the bubbles in each of the Figs. Illfa-c. and can 
therefore be obtained by cutting three different structures of orbit pairs. 

We have thus seen that by cutting through 2-encounters of the M(v [ -~' 2 ^ , 2) structures of orbit pairs considered, we 
obtain each family of <i-quadruplet once and once only, whereas each x-family is produced by three different structures. 
We therefore have 

, 2) = M d (v) + 3Af x (v) , (A9) 

and, by summing over v as in Subsection IA II 

-B m+ i =d m + 3x m . (A10) 
Eqs. QA8(1 and (|A10|) form a system of equations for the coefficients d m , x m , with the solution 

3 1 

d m = -A m + -B m+ i (Alia) 
x m = --A m - -B m+1 (Allb) 

3. Recursion relations 

We must calculate the auxiliary sums A mi B m defined in (|A1|) for both the unitary and the orthogonal cases. 
We start from the recursion for the number of structures Af(v) already used to evaluate the average conductance in 
Subsection III Dl 

M(v,2) -^A^ M ^ +1 U + 1) = (j-l)Af(^); (A12) 

k>2 ^" J 

see Eq. (O. This time, (|A"T2|) has to be multiplied not with {-l) v ^\ but with (-l) v ^L(v) = 
_(_l)n^' 2 ^ +1I )(£(#,2^fc+i]) + i) = -(-l)i / («' 2 ~ >1 )(L(t7[ 2 ^]) + 2). If we subsequently sum over all v with M(v) = m, 
our recursion turns into 

M(v)=m ( 

J2 \ (-l) v ^L(v)Af(v, 2) + £(_i)^ fc ' 2 -* fc+11 )(L(#.2-*+i]) + ij^fc.a-fc+i], k + 1) 

v [ k>2 

, M(v)—rn 

s l) £ (-l) V{v ^ ]) (L(v^) + 2W(v). (A13) 
Changing the summation variables as in Subsection III Dl we find 

M(v) = 

£ { (-i) v(vl L(vW(y, 2) + £(-i) y ^(L(«) + Wft k + 1) 

t? [ fe>2 
2 x M(v)=m-1 



i (~i) v{vl (m + w(v)- (Ai4) 



Using Af(v, 2) + £ fe > 2 A/"(v, k + 1) = Af(v), we can simplify the left-hand side to ( " )_m (-l) v(i7) (^(«) + WW ~ 
Y^rg^ - ™ 1 (— l) v ^J^f(v, 2) = A m — B m , whereas the right-hand side turns into — — lj (A m _i + c m _i) with c m _i = 
Y*? { -^ =m -' 1 (-l) v i^M(v), see Eq. We thus find a first relation between A m and B m , 

A m - B m = - - 1^ (A m _i + c m - x ) . (A15) 
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A second relation between A rn and B m follows from the recursion 12 

Af(v,3) - ]TAA(#^ fc+2 U + 2) = 2 (- - l) M(v^ 2 \2) + + W(v^) , (A16) 

k>2 \P J P 

which has to be multiplied with (-1) V W = _(_i)V(« [fe ' 3 ^ fc+2] ) = (_i)V(« [3 ^ 21 ) = _(_^v(v^) and summed over 
all v with M(v) = m. It is easy to see that the changed vectors yl k > 2 ^ k + 1 } } tJI 3 ^ 2 ^ t/I 3 ^] in the arguments of M have 
^■^fc,3-»fc+2]-j _ m ^ 7Vf(z/[ 3 ^ 2 ]) = 772 — 1 , and M(v^^) = m — 2. Again transforming the sums over v into sums over 
the arguments of N, we are led to 

M(v)=m ( 

£ (-l)^J^, 3) + ^^(77,^ + 2) 

v [ fc>2 

= 2 (^- 1 ) E (-1)^^,2)-- ^ (-l)^(L(fl) + 1)^(0). (A17) 

V V 

Eq. (T^TTt can be simplified if we use Af(v, 3) + X)fc>2 -^"(^ k + 2) = Af(v) — Af(v, 2), and recall the definitions of A m , 
B m and c m , to get 

cm - B m = 2 Q - 1^ B ro _i - |^ m _ 2 (A18) 
4. Results 

In the unitary case the two relations between A m and B m , Eqs. (|A15I) and l|A18|l yield a recursion for both 
quantities. In the unitary case, ^ — 1 = 0, c m = 0, the two equations simplify to 

A m - B m = 0, -B m = -A m _ 2 , (A19) 

i.e., A m and B m coincide and depend only on whether m is even or odd. Since in the unitary case d m — A mi x m — 
—B m+ i, see l|A3ll . and since the initial values d\ = 0, d% = 1 are already established, we come to the expressions (|4t)|) 

for Xm , d ra . 

In the orthogonal case § — 1 = 1, c m = (— l) m , Eqs. ljA15|) and (|A18fl take the form 

A m -B m = -A m _i + (-l) m (A20a) 
{-l) m -B m = 2S m _ x -2A m _ 2 . (A20b) 

Eliminating B m and -B m -i in (|A20b|) with the help of (|A20a(l . we find a recursion for A m , 

A m = — 3A m ^i . (A21) 

An initial condition is provided by A\ = —3, which accounts for the Richter/Sieber family of trajectory pairs with 
one 2-encounter, V = 1, L = 2 and m = 1. We thus obtain 



A m = (-3) m , B m An + A ro _! + = (-l) m (2 • 3" 1 - 1 - 1) . (A22) 

After that, using HAllaJI . I jAllbj l we arrive at the result 143(1 of the main text. 



Eq. IA16I follows from the special case I = 3 of Eqs. (42), (54) in [jj. To understand the equivalence to [jj, we need the iden- 
tity M(v^^\l) { = ] Af(^^) *P L( tf-Ti | 1 W (g' 3 ^ 1 ',l) = L( tZ] ]) J^(^ ] ) ( *= ) (i^ 3 "* 1 ) + l)^^ 3 -]), following from (*) 

A/'('0',1) = A^(?/'' 1— ) (see the footnote in Subsection 111 Pi . Eq. 1261 of the present paper, and (**) the fact that vi = and thus 

[3-*l] i 
v[ =1. 
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APPENDIX B: AN ALGORITHM FOR COUNTING TRAJECTORY PAIRS AND QUADRUPLETS 

In our approach transport properties phenomena are related to trajectory pairs or quadruplets differing in encoun- 
ters; we showed that all topologically different families of pairs or quadruplets can be found using the corresponding 
structures of orbit pairs considered in the theory of spectral fluctuations The number of structures to be consid- 
ered grows exponentially with the order of approximation, and in absence of an analytic formula (in particular, for 
Ericson fluctuations or crossover in a magnetic field) can be evaluated only with the help of a computer. Here we 
describe an algorithm for the systematic generation of orbit pairs. 

Assume time-reversal invariant dynamics, and let A be a periodic orbit and B its partner obtained by reconnections 
in a set of encounters associated to the vector v — («2, ^3, ^4, • ■ ■)■ us then number all the L = L(v) encounter 
stretches in A in order of traversal. Each encounter has two sides which can be arbitrarily named left and right such 
that the two ends of the i -th encounter stretch can be called its left (it) and right (i r ) ports. The direction of traversal 
of the encounter stretches in A will be denoted by another vector ~~o = {o\Oi . . . ctl). Here Oi is equal to 1 if the i -th 
stretch is traversed from the left to the right such that ii and i r are its entrance and exit ports respectively, and equal 
to -1 in the opposite case. There can be 2 L different a. However, only relative direction of motion within encounters 
is physically meaningful: we can, e.g., assume that the first stretch of each of the V (v) encounters is passed from 
the left to the right. This leaves 2 L ~ V physically different a. Each orbit link connects the exit port of an encounter 
stretch of A with the entrance port of the following stretch. Both of these ports can be left or right which makes four 
combinations possible, i r — ► (i + i r — ► (i + l) r , ii — > (i + l)i, ii — > (i + 1);, the sense of the link traversal in 

A being in all cases i — ► i + 1. These choices are uniquely fixed by the vector a. 

In the partner orbit B the left port ii will be connected by an encounter stretch not with the right port i r but with 
some other right port f(i) r - The set of reconnections of all ports can be written as a permutation 



in which the upper and lower lines refer to the left and right encounter ports, correspondingly. Since reconnections 
are possible only within encounters, P e nc must consist of as many independent cyclic permutations (cycles) as there 
are encounters in the orbit A, with each cycle of / elements corresponding to an Z-encounter. The links of the orbit 
B are unchanged compared to A, but may be passed with the opposite sense. The orbit B may exist in two time 
reversed versions; we shall choose the one in which the encounter stretch with the port 1/ at its left is passed from the 
left to the right. This choice made, the sequence of visits of all ports and the direction of traversal of all encounter 
stretches and links in B become uniquely fixed by the permutation P cnc and the the vector ex. Indeed, let us start 
from 1;, move along the encounter stretch arriving at the right port f(l) r and then traverse the link attached. Where 
we move next in B depends on the ports connected by this link in the original orbit A: 

1. /(l) r — > [/(I) + 1]/ ■ The sense of traversal of the link in B is the same as in A. The next encounter stretch is 
traversed from the left to the right leading from the port p' 2 = [/(l) + l]i to [f(p 2 )] r - 

2. /(l) r — * [/(l) + l] r ■ The sense of traversal of the link in B is the same as in A. The next encounter stretch is 
traversed from the right to the left leading from p' 2 — [/(l) + l] r to [/^(Pa)],; i-e., to the element of the upper 
row in P cnc corresponding to f(l) + 1 in the lower row. 

3. [/(I) — 1] ; — > /(l)r • The link is traversed in B in the direction opposite to A leading from /(l) r to p' 2 = 
[/(l) — The next encounter stretch leads from the left port p' 2 to the right port [/ (pQL ■ 

4. [/(I) — l] r — > /(l)r- The link is traversed in B in the direction opposite to A leading from /(l) r to p' 2 = 
[/(l) — l] r . The next encounter stretch leads from the right port p 2 to the left port [Z -1 ^)]/- 

Continuing our way we eventually return to the starting port 1/. If that return occurs before all 2L encounter ports 
are visited the combination P onc , a has to be discarded since it leads to a partner consisting of several disjoint orbits, 
a so called pseudo-orbit. Otherwise we have found a structure of the periodic orbit pair (A, B) with the encounter 
set v and established the port sequence in B as well as the sense of traversal of all its links and encounter stretches. 
Running through all LA/ JIz>2 l" lv ^ permutations associated to v and through all 2 L ~ V essentially different a we find 
all N(v) structures of the orbit pairs. 

Suitable cuts yield all trajectory doublets and quadruplets relevant for quantum transport. E.g., if we cut A and B 
in the initial link (i.e. in the link preceding the port 1;) we obtain a trajectory pair contributing to the conductance; 
the numbers jii , p a needed for calculating conductance in a magnetic field are obtained by counting the number of 
links and encounter stretches changing their direction in B compared with A. 




(Bl) 
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The trajectory quadruplets contributing to the conductance covariance and other properties are obtained, in ac- 
cordance with our method I (see Appendix IXI) . by cutting the pair (A, B) twice, in the initial link and in any of the 
L orbit links producing thus L + 1 quadruplets per orbit pair. If the second cut goes through a link preserving its 
sense of traversal the result is a <i-quadruplet, otherwise it is an ^-quadruplet. (We remind the reader that in the last 
case the sense of traversal of the trajectories 7, 8 is changed to the opposite compared with the periodic orbits; this 
has to be taken into account in calculation of the transport properties in the magnetic field.) Method II of producing 
quadruplets consists of cutting a 2-encounter out of the orbit pair (A, B) which is possible only if the encounter set 
v contains at least one 2-encounter, i.e. «2 > 0; the resulting quadruplet will be characterized by the encounter set 

The relatively trivial case when time reversal is absent can be treated by choosing <jj = +1, i = 1,...,L (all 
encounter stretches are traversed from the left to the right, and all links are attached to the ports like i r — > (i + 1) ; ). 




APPENDIX C: DIAGRAMMATIC RULES FOR ARBITRARY MULTIPLETS OF TRAJECTORIES 

Our semiclassical techniques can be expanded to a huge class of transport problems, for cavities with arbitrary 
numbers of leads, and observables involving arbitrary powers of transition matrices. We then have to evaluate (sums 
of) general products of the type 

\m— 1 / 

with a± , . . . , <Zj \ , ■ ■ ■ , ^ denoting 2M mutually different channel indices associated to any of the attached 
leads, and II a permutation of 1, 2, ... , M. Using the semiclassical transition amplitudes JQ), we express (Z) as a sum 
over multiplets of trajectories a.k,/3k connecting channels as a m (aj m ' ) — > a 2 " l ' ) ), Pmi^i^ —> )j 

^^...^6^-+-+^-^--^)/^ . (C2) 

Contributions to (|C2|) arise from multiplets of trajectories where f3i, . . . , (3m either coincide with ai, . . . , olm-, or differ 
from the latter trajectories only inside close encounters in phase space. (Z) thus turns into a sum over families of 
multiplets characterized by a vector v. 

Proceeding as in Subsections III CI and IIII Bl we represent the contribution of each family as a sum over the trajec- 
tories ax, ... , ctM and an integral over the density of stable and unstable coordinates, 

(Z)\ fam = ^( £ |A^.-KJ 2 /^ s ^^^ ; 

H \ai,...,a M ' / 

(C3) 

here, w(s, u) is obtained by integrating {tt L ~ v Yi a =i ^cnc( s i over the durations of all links except the final link of 
each of the M trajectories. Doing the sum over ctk with the Richter/Sieber rule, we find the same link and encounter 
integrals as before, and thus a contributions from each link and a contribution —N from each encounter. 

The same rule applies for products slightly different from ijClfl . First, if some of the channels a[ , . . . , a^ M ^ or some 
of the channels , ■ ■ ■ , a>2^ coincide (as in many examples studied in the main part) some of the subscripts in (|C1|) 
will appear not twice but 4, 6, 8, . . . times. We then have to consider all possible ways to pair these subscripts. In the 
spirit of Wick's theorem, each of these possibilities contributes separately. (Note that, if a subscript appears an odd 
number of times, the corresponding product cannot be related to multiplets of trajectories, and may be expected to 
vanish after averaging over the energy). Second, in the orthogonal case the two subscripts of one t or t* in IjClJI may 
be interchanged without affecting the final result; the corresponding trajectory is then reverted in time. 

With these rules, one can evaluate a huge class of observables relevant for quantum transport. For each single 
application, only the counting of families remains to be mastered. 



APPENDIX D: SPECTRAL STATISTICS REVISITED 



We here want to reformulate our previous results on spectral statistics [9| in the present language of diagrammatic 
rules. In contrast to @, we start from the level staircase N(E), defined as the number of energy eigenvalues below E. 
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N(E) can be split into a smooth local average N(E) and an fluctuating part N osc (E) describing fluctuations around 
that average. We want to study the two-point correlation function of N osc (E) 

^) = (n osc {e + ^)n osc (e-^)). (Dl) 

The latter correlator yields the spectral form factor K{t) — i dee 2leT R(e) through the identity R(e) = 
ij {^^p [E + 2^ {e - 2^pj"j = -n 2 ^-- To check on the last member of the foregoing chain of equa- 

tions one must write the average (•) in (|D1(I as an integral over the center energy E, take its second derivative by e 
and integrate by parts the terms containing N"N in the integrand. 

In the semiclassical limit, Gutzwiller's trace formula determines N osc (E) as a sum over periodic orbits A of arbitrary 
period T A (E), N osc (E) = — Im^^i^e 1 ^^; here, F A depends on the stability matrix M A and the Maslov index 
\i A of A as F A = 1 e 1 ^- 4 % , and SAE) is the classical action of A at energy E. The correlation function 

y/\ det(AfA-l)| 

C(e) turns into a double sum over periodic orbits A and B, 

C( £ ) = -^Re (j2 F A F^ s ^- s ^y h e i TM % TBlE) e ) , (D2) 



\A,B 



where we have used S A (E + « ^(E) + ^(E)^, S 8 (£; - « S B (£) - T B {E)^= and T H = 2nhp. 

To evaluate the contribution to l|D2l) resulting from a given structure of orbit pairs differing in encounters (see 
Subsection III ( J[) . we replace the sum over B by an integral over a density w(s, u) of phase-space separations inside A. 
Similarly as for transport w(s,u) is defined as the integral of riL _ vr ,v 1 ta , — 7 over all piercing times; integration 

over the first piercing time leads to multiplication with the orbit period, whereas the remaining integrals can be 
transformed into integrals over all link durations but one. Approximating Fg « F A and T& ss T A , we find 

= a^I (E l^l 2 / ^ s / ^-^^(«^)e i(S - (E) - SB(E))/R e 2ir ^/ r ^ . (D3) 

We here divided out L, because for each orbit pair any of the L links may be chosen as the "first" ; without this 
division, each pair would be counted L times. The sum over A can now be done using the sum rule of Hannay and 
Ozorio de Almeida [45j |. 

J2\F A \ 2 S(T-T A ) = ^ => ]T|F4| 2 (.)= jdrU-). (D4) 

A A J 

The multiplication with the orbit period is thus replaced by integration over the period or, equivalently, over the 

duration of the remaining link. Writing T A — J2f =1 U + Y^=i ^enc( s i M )> we can now s P ut JESt mto * ne prefactor 
„ \ T , an integral 

27T- L 1 ° 

/ dUe 2 ^" = -±£- (D5) 



for each link and an integral 



d'-W- 1 ! ^,,^ - - e a'.^M^H e i£ i ..^/t) = £££ (D6 ) 



2Ue 

"fi 1 - 1 ' / ~ 7^ 

for each encounter; to ensure convergence we assume e to have an infinitesimal positive imaginary part. Since all Tu 's 
cancel, we obtain a factor —5^ for each link, and a factor 2Zie for each /-encounter. The overall product reads 

^)L m = ^(-^ L W vlR <^ v - L - (° 7 ) 
1 

The corresponding contribution to the spectral form factor is easily evaluated as 

Vf^\\ (~l) y IL^' L-V+l ;rvo\ 

K{T) \^ = {L~V-l)\L T ■ (D8) 
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We have thus rederived one of the main results of [9| in the elegant fashion suggested by the present work on transport. 
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